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ANALYSIS OF THE INJECTION OF A HEATED, TURBULENT JET INTO A MOVING 
MAINSTREAM, WITH EMPHASIS ON A THERMAL DISCHARGE IN A WATERWAY 

by 

J ames F. Campbell 
(ABSTRACT) 

An experimental, and. theoretical, investigation has been undertaken 
to study the trajectory and growth of thermal effluents having a range 
of discharge velocities and temperatures. The discharge of a turbulent 
effluent into a waterway was mathematically modeled as a submerged jet 
injection process by using an integral method which accounts for 
natural fluid mechanisms such as turbulence, entrainment, buoyancy, and 
heat transfer. The analytical results are supported by experimental 
data and demonstrate the usefulness of the theory for estimating the 
location and size of the effluent with respect to the discharge point. 
The capability of predicting jet flow properties, as well as two- and 
three-dimensional jet paths, was enhanced by obtaining the jet cross- 
sectional area during the solution of the conservation equations (a 
number of previous studies assume a specific growth for the area). 
Realistic estimates of temperature in the effluent were acquired by 
accounting for heat losses in the jet flow due to forced convection and 
to entrainment of free-stream fluid into the jet. 
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VI . INTRODUCTION 


One of the most basic processes in fluid mechanics is the mixing 
that occurs during the interaction of two intersecting stre ams of 
fluid. A large portion of interaction processes exists where one 
stream of fluid has a much smaller mass (or volume) flow than the 
other stream, such as the case of a plume issuing from a smoke stack 
or a fluid injecting into a boundary layer on a vehicle traveling at 
hypersonic speeds. This class of interaction problems can be cate- 
gorized as fluid, or jet, injection processes, and are usually asso- 
ciated with a nearby solid boundary or surface, and are complicated 
by the fact that in their natural state they are almost invariably 
turbulent . 

An example of a fluid injection process of current interest is 
where a heated water effluent is discharged into a river, reservoir, 
or estuary , the water being initially heated during the course of 
some industrial process. For example, the generation of electricity 
accounts for approximately TO percent of the waste heat discharged 
daily into waterways in the nation.^ In order to satisfy the nations' 
demand for electricity in 1980, almost 250 billion gallons of water per 
day will be required for cooling, of which about 200 billion gallons 
will be from fresh water sources. This will comprise 1/5 of all fresh 
water run-off in the country. A schematic diagram is presented in 
Fig. 1 (Aronson ) to illustrate how fuel in an electric power plant 
becomes electrical energy and how water from a nearby river or reservoir 


1 












3 


becomes heated. The steam used to turn the turbines and hence to 
generate electricity must be condensed in order to be re-used. Cooling 
water drawn from a nearby water source circulates through the condenser 
withdrawing heat from the steam, and in its warmed-up state is then 
returned to the water source. The increase in cooling water tempera- 
ture above the ambient temperature is typically 10° to 25°F, depending 
of course on the volume flow of water through the condenser and the 
amount of heat extracted. This process of heat disposal, called the 
"once-through" cooling method, is the most widely used technique 

because it is the easiest to build and the cheapest to operate 

2 

(Woodson ). At present about 90 to 95 percent of steam-electric cooling 
water is put back into the waterway using this once-through method; 
the remaining percentage of water is treated in cooling towers and/or 
ponds before discharge. 

In an attempt to control the growth of this problem, water 
quality standards have been legislated which set temperature limits 
(above the ambient temperature) on the various types of water bodies 
accepting heated water. All of these standards apply to that region 
of the receiving water outside of the "mixing zone", where the heated 
effluent mixes with the ambient water. The fact that the definition 
of the mixing zone varies not only from state to state, but often 
between different locations within the same state, adds to the ambiguous 
nature of the problem. 

A great deal of information has been generated concerning the 
effects of heat on aquatic life (Refs. 4 to 6), however it has only 
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been in recent. years that serious research attention has been turned to 
the fluid physics involved with modeling the heated discharge. Many 
of the previous studies directed towards the physical understanding of 
the calefaction of waterways have been more concerned with the far-field 
effects of the thermal discharge. Examples of this are the works by 

7 

Edinger and Geyer who consider the overall heat balance of bodies of 

Q 

water having excess heat addition, and Mahgary who attempts to examine 
the diffusion of water temperature assuming a two-dimensional mixing 
model, but neglects to consider momentum transport. Since it is 
desirable to observe the time history of a heated effluent (or thermal 
plume) from the point of discharge to the point of final thermal, mixing 
with the ambient fluid, it is therefore necessary to consider the near- 
field effects. In their general reviews of research needs, Tichenor^ 
and Parker 10 stressed the importance of considering the momentum jet as 
the method of modeling the discharge in the near field. This is 
certainly the more pragmatic approach to the physical understanding of 
the problem, particularly if the definition of a thermal mixing zone is 
desired in order to minimize impact on the local water environment. 

Accordingly, the present study was initiated to investigate the 
pertinent fluid mechanics and heat-transfer aspects of heated, turbulent 
effluents discharging into a moving water environment. Emphasis has 
been given to determining the effects of the ambient flow variables and 
effluent discharge properties on the three-dimensional path, or tra- 
jectory, of the thermal plume and to examining how the flow properties 
of and heat loss from the plume vary along the path. In a realistic 
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situation, the heated effluent would he discharged into an irregular 
velocity field typical of most river flows, however, in order to pro- 
vide a starting point for the present theoretical analysis, the 
effluent is considered as a submerged jet flow injecting into a moving 
water system having a spatially dependent velocity field. The general 
development of the theoretical model is contained in Chapter VIII, with 
details contained in the Appendices in Chapters XIV and XV. The theory 
is compared with other analytical models in Chapter X, as well as with 
experimental data acquired from a number of investigations. Experimen- 
tal tests were conducted for the lateral, vertical, and oblique jet 
injection into a water channel, the Jet having a range of injection 
velocities and temperatures . Results of the tests are described in 
Chapters IX andXIII and axe presented to verify the major assumptions 
of the mathematical model. Results are represented by properly non- 
dimensionalized parameters to make the analysis applicable to any 


situation. 



VII. REVIEW OF LITERATURE 


oince the question of siting thermal discharges on a moving water- 
way will largely be settled on the basis of localized effects near the 
outfall, it is important for us to understand the fluid mechanics of the 
jet injection process. Margason's bibliography 11 on previous experi- 
mental investigations of incompressible jets injecting into a cross 
flow provides a starting point for the discussion. The complex inter- 
action that takes place after a fluid is injected into a moving stream 
results in a three-dimensional flow field, even though the Jet may be 
following a two-dimensional path. This, of course, complicates the job 

of the experimentalists. Gordier , Margason , and Platten and 

lh 

Keffer , concentrated their efforts on measuring the jet's trajectory 
for a wide range of injection angles (a.) and velocities (V./V ), the 

1 1 CO * 

trajectory being one of the easiest jet properties to measure (see 

Fig. 2). Jordinson ^ , Keffer and Baines 1 ®, Ramsey 1 "^, and Kamotani and 

l8 

Greber measured details of the interaction process and revealed the 
very complex flow field that exists due to the pressures and shear 
stresses in and around the jet flow. Very good descriptions of how 
the jet flow distorts and develops under the influence of the free- 

stream flow and body forces are given by Keffer and Baines 1 ®, Ramsey 1 ^, 

19 20 

Abramovich , and Keffer 

The most detailed measurements of the flow field appear to be those 
of Keffer and Baines 1 ®, and Kamotani and Greber 1 ®. In particular, 
Kamotani and Greber examine the structure of the rotational velocity 


6 
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field which results because of the shearing action between the free- 

stream flow and the edge of the jet flow. The rotational velocity 

field, usually interpreted to be a pair of counter-rotating vortices, 

is an interesting as well as important facet of the injection process; 

it affects the path of the jet as well as the mechanisms that govern 

entrainment. Both Refs. 16 and 18 reduced their data so a measure of 

mass flux in the jet could be obtained, which in turn resulted in 

estimates of the rate that free-stream fluid is entrained into the jet 

structure. In addition to investigating the characteristics of the 

IT 18 

velocity field, Ramsey and Kamotani and Greber measured the 

temperature field resulting from the injection of a heated Jet. 

It is important to note the limiting conditions of a jet injected 

at an arbitrary angle into a cross flow. If the angle between the jet 

axis and the direction of the free-stream velocity (a) goes to zero for 

a given V^/V^, the condition exists of a jet in a co-flowing stream. 

If the jet velocity becomes very large (V^/V^ ■+■ 00 ) for a given a, the 

situation approaches that of a free jet. An indication of the magnitude 

of entrainment for a jet in a co-flowing stream and for a free jet has 

21 22 

been provided by Morton , and Ricou & Spalding , respectively. 

In reviewing some of the theoretical methods available for modeling 
an injection process, it should be emphasized that a method is desired 
here which allows the thermal effluent to be followed from its point 
of discharge to some point downstream where complete thermal mixing 
takes place. At a minimum, the method should be able to estimate the 
three-dimensional (3-D) path of the thermal plume, show how the plume 
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size varies as it procedes downstream, and include basic heat-transfer 
mechanisms which allow the plume temperature to be predicted. The most 
advanced method would result in a complete and detailed description of 
the flow field resulting from an injection into a moving stream. This 
could be done by solving the full three-dimensional , turbulent , Navier 
Stokes equations in an Eulerian framework, which requires the specifica- 
tion of the eddy viscosity field. 

In order to avoid the complexities inherent with this approach, 
many studies have tried to theoretically model the gross features of 
the injection process by describing the fluid motion of the Jet from 
the point of discharge in a Lagrangian framework. This procedure allows 
an estimate of Jet properties to be obtained if the appropriate forces 
acting on the Jet flow are accounted for. Since the jet path is the most 
obvious of the jet properties , it is natural that early attempts were 

19 

concerned only with obtaining estimates of the trajectory. Abramovich , 
for example, obtained the trajectory of a jet which had a circular 
cross-section at the injection point by balancing the centrifugal 
and blockage forces perpendicular to the trajectory. His basic argument 
was that the blockage effect of the jet flow on the free-stream flow 
could be approximated by assuming that the jet flow acts as a "solid" 
body inclined at some angle to the free stream. He accounted for the 
deformation of the jet's cross-section by assuming the shape to be ellipti- 
cal and by specifying a growth rate for the cross-sectional area, which 
is necessary if only one force equation is used to obtain a solution 
for the trajectory. One of the serious drawbacks of this method is the 
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assumption that the component of jet momentum perpendicular to the 
direction of the free-stream flow remains constant. This assumption 
was relaxed by Schetz and Billig . 

Two other forces acting perpendicular to the jet trajectory and 
which help govern the development of the jet flow are: (l) buoyancy 

force resulting from a difference in density between the Jet and free- 
stream fluids, and (2) entrainment force resulting because of the free- 
stream fluid that is drawn into the jet structure. Theoretical 
trajectories were obtained by Reilly 2 ^ and Campbell and Schetz 25 using 
procedures similar to those described above, but also accounting for 
the entrainment phenomenon; Campbell and Schetz, in addition, included 
the buoyancy force in their model. 

Since all of these previous works utilized an assumed area growth 

based on experimental data obtained in the proximity of the injection 

point, they are not suitable for providing realistic trajectory 

information further downstream. This deficiency can be avoided if, 

instead of assuming an area growth rate, a momentum conservation 

equation in the direction of the jet path is used; this equation is used 

in addition to the conservation equation perpendicular to the trajectory. 
26 27 28 

Wooler et. al. and Hoult et. al. * used this procedure, solving 
force equations normal and parallel to the jet path simultaneously 
to obtain a solution for the trajectory. An added advantage of using 
these two force equations is that, if all of the appropriate forces are 
accounted for, the solution procedure allows the jet flow properties to 
be estimated. The theoretical studies discussed up to this point, 
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however, do not include all cf the necessary forces. This is illus- 
trated in Table I which presents a summary of the capabilities of a 
number of theoretical investigations. In particular, the shear stresses 
that exist in the flow field due to differences between the jet and 
free-stream velocities are not represented in most theoretical efforts. 
The shear force normal to the trajectory is usually combined with the 
force due to the pressure distribution around the jet flow to form the 
blockage force (mentioned previously). The shear force parallel to the 

jet path must be accounted for in a force balance on the jet flow, 

29 30 

which was done by Hirst and Campbell and Schetz 

The theories discussed thus far have been concerned with predict- 
ing the flow characteristics of a jet following a two-dimensional (2-D) 
path. This type of path occurs when the radius-of-curvature vectors 
associated with the trajectory lie in one plane, the trajectory of 
course being a curve in that plane. If the injection and free-stream 
velocity vectors shown in Fig. 2 are considered to form an "injection 
plane", then a two-dimensional trajectory will always result when this 
plane is oriented vertically, i.e. aligned with the gravity vector. 

When the injection plane is rolled away from the vertical, a three- 
dimensional trajectory may occur depending on the buoyancy of the jet 
flow. If the buoyancy force is absent, or is small with respect to the 
jet momentum, then a two-dimensional path will result. A larger 
buoyancy force, however, will cause the jet to bend out of the injec- 
tion plane, and the result will be a three-dimensional trajectory. 

The radius-of-curvature vectors associated with this type of trajectory 



TABLE I. - CAPABILITIES OF PREVIOUS THEORIES USING INTEGRAL TECHNIQUES 



Forces parallel and perpendicular 
to jet trajectory 


Investigators 

Arbitrary 

Injection 

Arbitrary 

Injection 

Blockage 

Entrainment 

Buoyancy 

Viscous 

Shear 

Jet Flow 
Properties 

3-D 

Trajectory 


Velocity 

Angle 





P .A,V 

T 


Abramovich^ 

X 

X 

X 

■ 
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Schetz and Billig 

X 

X 

X 
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Reilly 

X 


X 

X 
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Campbell and Schetz 

X 


X 

X 

X 





Wooler et . al . ^ 

X 


X 

X 






27 

Hoult and Fay 

X 



X 

X 





Pfi 

Hoult and Weil 

X 



X 

X 





29 

Hirst * 

X 

X 

X 

X 

X 

X 

X 


X* 

30 

Campbell and Schetz 

X 

X 

X 

X 

X 

X 

X 

X 


Present 

X 

. X 

X 

X 

X 

X 

X 

X 

X 


# No 3-D trajectory results given 
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do not lie in one plane. Cince u Jet with a two-dimensional trajectory 
is a special case of the Jet with the more general three-dimensional 
trajectory, it is desirable to have a theory that can estimate the more 
general situation. At present, the only theoretical method which 
appears capable of predicting jet flows with three-dimensional tra- 
jectories is that of Hirst, although no results of this kind were 
presented in his report (Ref. 29). 

In view of the preceding comments , the purpose of the theoretical 
portion of the present investigation is to develop an integral method 
which accounts for natural fluid mechanisms such as turbulence, 
entrainment, buoyancy, and heat transfer in the conservation equations 
governing the jet flow. In particular, it is desirable to have a 
theory that: (l) utilizes the momentum conservation equation along the 

Jet path, in addition to that perpendicular to the path, in order to 
avoid any assumption regarding the growth of the jet's cross-sectional 
area as was done in Refs. 19, 23-25; (2) accounts for a shear force in 
the momentum equation along the jet trajectory, which was not included 
in the works in Refs. 26-28, in order to provide improved predictions 
of two-dimensional jet trajectories, and to allow estimates of jet 
flow properties, such as velocity, cross-sectional area, and momentum, 
to be made; (3) obtains a third momentum conservation equation for the 
jet flow which, when solved simultaneously with the other two momentum 
equations, allows three-dimensional jet paths to be calculated, thus 
extending the theories reported in Refs. 29 and 30; (4) provides 
estimates of the temperature of the Jet fluid by examining several heat 
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transfer mechanisms that can account for the heat loss from the jet 
flow; and (5) is easily adapted to account for free-stream flows with 
either a nonuniform velocity field or a nonuniform temperature field. 
The purpose of the experimental portion of the present study is to 
provide information for a submerged jet injected laterally, vertically, 
and obliquely into a water channel in order to support the theoretical 
effort , as well as to extend the experimental results of previous 
investigations. It is desirable to examine the effects of changing 
injection velocity for the different injection orientations, and to 
provide test results at higher injection velocities than those reported 
in Refs. 12 to 20. In addition, the effects of increasing injectant 
temperature on the injection processes should be ascertained, and 
measurements of jet temperature decay made to compare with the results 
of Refs. 17 and 18. 



VIII. THEORETICAL DEVELOPMENT 


This chapter is concerned with the development of a theory which 

approximates the fluid mechanical process that occurs when a submerged, 

turbulent jet of circular cross-section is injected into an effectively 

infinite free-stream flow. The mathematical model should allow the 

jet to penetrate into the cross flow , to bend over and spread under 

the influence of natural fluidic forces, and the jet's velocity vector 

to approach the free-stream velocity vector at some point downstream 

from the jet exit. This capability is obtained by considering a 

section of jet fluid as a control volume similar to the approach used 
2b 

by Reilly . This is illustrated in Figure 2 which depicts the 
trajectory of a jet injecting into a free-stream flow having a velocity 
V taken to be spatially but not time dependent. The origin of the 
cartesian (x,y,z) axis system is at the injection point, while the 
natural (s,n,t) axis system moves and rotates as it follows the path 
of the jet axis which is traced out by the jet velocity V. The s-axis 
is located along the trajectory, while the n-axis is oriented perpen- 
dicular to the trajectory in the direction of the radius of curvature 
of the trajectory; the t-axis is perpendicular to both the s- and n- 
axes . A complete discussion of the natural coordinate system is given 
in Appendix B. The equations expressing conservation of mass and 
momentum in the direction of the natural axes are derived in the 
following sections. 
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Figure 2.- Cartesian and natural coordinate systems. 
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Continuity 

The following integral expression equates the net influx of mass 
into the control volume to the rate of increase in mass in the control 
volume : 


it JfJ p£ da = -f f ps7 * V 6 (1) 

where d6 and da represent the elemental surface area and volume, 
respectively. It is assumed that the flow process is steady (in the 
mean), fully turbulent, and incompressible. The fact that we assume 
incompressibility does not imply that the flow process is one of 
constant density. Equation (l) thus becomes, 



Carrying out the operations suggested by this equation leads to the 
mass flows through the surfaces .of the control volume and 

© in Fig. 2). 

It is noted that d6 = dA for surfaces © and (2), which 
represent a cross-section of the jet flow; i.e. areas that are 
perpendicular to the trajectory. Since mass flow is a continuous, 
single-valued function of position along the trajectory, a Taylor 
expansion can be performed to obtain the mass flow through surface 
as a function of the mass flow through surface (T) . The difference 
between mass flows through and (2) represents the mass flow through 
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the sloping face ©■ which defines the amount of free-stream fluid 
mass (per unit of jet length) that is entrained (E), or drawn, into 
the control volume. This can be written as. 



The two most common profiles used in the literature to describe 

the velocity variation at a given cross-section in the Jet flow are 

the Gaussian and Top-hat. The Gaussian representation is particularly 

useful for providing estimates of flow properties on the Jet's 

29 

centerline and was used in the development of Hirst . In using the 
Gaussian profiles, however, the assumption is made that the jet flow 
is circular, i.e. axisymmetric , which is valid for free-Jet flows and 
for co-flowing jet flows, but is of questionable use for a jet 
injecting into a cross flow, where the jet cross-section is not 
circular"^. In addition, the theoretical results obtained by using a 
Gaussian velocity profile are only applicable in the region where the 
jet flow has become fully developed. 

Top-hat profiles represent the average jet flow properties and 
have been used in a number of theoretical studies * * * . By 

using the average flow properties in the conservation equations, it is 
not necessary to place restrictions on the symmetry of, or to assume 
similarity in, the jet flow. This means that the governing equations 
can be used to describe the jet flow in the region where the jet is 
fully developed, or in the region where the jet is only partially 



18 


developed and a potential core still exists. Accordingly, the present 
study uses the averaged jet flow properties defined as follows: 



which state that at a given location along the trajectory the jet's 
local velocity, density, and temper&ture values are integrated over 
the jet cross-sectional area. This permits the conservation of mass 
in the jet, equation (3), to he expressed in differential form as, 

21 2k 2& 

Morton , Reilly , and Wooler et. al. obtained continuity equations 
of this form. 

In recent years it has been recognized that inclusion of the 
entrainment process in a jet injection analysis is important not only 
because of its influence on jet momentum, and hence trajectory, but 
also because this process allows for the mixing of jet and free-stream 
scalar properties such as temperature and salinity. The fact that 
entrainment occurs when there is relative motion between two flow 
fields has been used by different researchers to justify relating 
entrainment to the appropriate velocities normal and parallel to the 
jet axis. This type of representation of the entrainment function (E) 
results in a variety of empirical "constants" which must be adjusted 
in order to obtain suitable agreement with experimental data. Examples 
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of several attempts to approximate E in this fashion are seen in the 

works of Kamotani & Greber 18 , Wooler et. al. 28 , and Hoult et.al. 2 ^* 28 

The model of entrainment presented by Wooler et. al. uses three 

constants, one of which is obtained by satisfying the Ricou-Spalding 
22 

measurements for an isothermal free-jet. Experimental trajectory 

data, rather than mass flux information, were used as the criterion 

for indirectly adjusting the other two constants. Experimental mass 

flow data are presented by Kamotani & Greber to support the constants 

used in their model of the entrainment function. 

Because of the complex helical streamline pattern (us ual ly 

interpreted as a pair of counter-rotating vortices) evident in the 

lee side of the jet, it is believed that the entrainment function 

cannot be split into two totally independent parts as assumed in the 

above description. One attempt to account for the free-stream fluid 

entrained into the jet as a result of the helical circulation pattern 

31 

has been reported by Platten & Keffer . Their entrainment model was 
29 

extended by Hirst to account for local buoyancy in the jet flow. 

In consideration of the complicated nature of analytically 
predicting the entrainment function and the desire to keep empirical 
constants to a minimum, the present study defines this jet property 
by using the experimental data for an air jet obtained by Keffer & 
Baines^ 8 . Equation (6) represents 

E * § PJS»<V - VJ 


( 6 ) 
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the entrainment parameter in functional form, where the entrainment 
coefficient (E*) was obtained from measurements of mass flux along the 
4et axis 1 ^. In order to use this formulation for E, V must always be 
equal to or greater than V^. The entrainment coefficient is presented 
in Figure 3 for three injection velocities and, as noted, is a 
function of both injection velocity and distance along the trajectory. 
An empirical expression was obtained here to represent E* in the 
present mathematical model and is shown compared with the experimental 
data. It is assumed that this function can be used at larger values 
of VR and s/d^ than those shown in the figure. 

It should be noted that the empirical expression for E® implies 
that as s -*■ 0, E* (and hence E) -► 0. This is an unsettling 
possibility because, near s » 0, the Jet flow in this type of injec- 
tion process might be expected to resemble that of a free jet, having 
similar entrainment rates. Keffer S> Baines 1 ^ observed that as s 0 
E* was of the same order as that found by Morton 21 for simple jet flows. 

This is further substantiated if we consider the work of Ricou & 

22 

Spalding who measured entrainment for axisymmetric turbulent free- 
jets. The formula they suggested as best representing their data is, 


pAV 

T^77 


0.32 


d. 

x 


(7) 


This can be differentiated with respect to s to get the entrainment 
in the free jet. Putting the results in the same form used by 
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Experiment (Keffer and Baines^) 



Figure 3.- Entrainment coefficient as a function of position 
along axis of isothermal air jet; <x = 90°. 
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Morton we get , 


E = 0.08 C.p.V. (8) 

ill 

Thus the constant 0.08 compares favorably with Morton’s value of 0.1l6. 

It is assumed in the present study that E* = 0.08 until the 
empirical expression presented in Figure 3 predicts a greater value, 
which is then used. This is reasonable in view of the fact that an 
increase in VR decreases E# so that as VR », E* should approach 
the free- jet value. 

Conservation of Momentum 

The following integral expression equates the rate of increase of 
momentum in the control volume to the sum of forces acting on the 
control volume plus the net influx of momentum into the control volume: 

it JJfv do - fff\ da d6 -/JW • t. d{) 

(9) 


where the surface stress tensor (x) is taken to contain both shear 



and pressure terms and f^ = body force per unit volume. Since we 
have assumed a time independent flow process , the term on the left- 
hand side of the equation is zero; hence. 




'N 


d<5 - 



( 10 ) 
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where the equation will be shown to yield expressions for momentum in 
the direction of the natural coordinates. 

Assuming the body forces are independent of the integral over do, 
we write 




( 11 ) 


where, as a first approximation, the volume is assumed to be equal to 
Ads, which neglects the rate of change of A with s. Since the 
components of f^, and later V^, are desired in the s, n, and t 
directions, we can use a vector identify to obtain. 


f = ( f 
x b u b 


e s )e s + (f b 


-r x-r ,-r 

e )e + (f 
n n b 


't )e t 


( 12 ) 


and 


~ /T 

V * (V 

oo a 


e )e + (V 

S3 « 


e )e + (V 
n n « 


e t )e t 


( 13 ) 


where V = V e . 

oo oo x 

There are two body forces considered to be acting on the jet 
flow. The first force is due to the buoyant condition which results 
from a difference in density between the jet and free stream fluids. 
This force acts in the y direction and is given by. 


F v = F, e = g(p m - p)e 
b by 00 y 


(lh) 
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The second force results from centrifugal effects associated with the 
jet having mass and following a curved path. This force acts in the 
direction of the radius of curvature of the trajectory and is 
expressed aB , 


F 

c 


F e 
c n 


m V 2 -* 
R e n 


(15) 


The total hody force on the control volume (f^) is the vector sum of 
the buoyancy and centrifugal forces and is substituted into Eq. (12) 
to obtain. 


F, (e 
b y 


e o )e + [F (e • e ) + F ]e + F, (e 
s' s b v y n c n b v y 


5 J 0 

v t 


(16) 


The pressure portion of the shear-stress tensor term in Eq. (10) 
can be written in component form as , 



(IT) 


where p^ represents the local pressure on the respective surfaces 
indicated by the integration. 


n- Momentum 


The n-momentum equation is obtained by taking the n components 
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of the various vector quantities in the momentum integral equation 
(10). The term representing the body forces, for example, is extract- 
ed from Eq. (l 6 ) as. 


f b n = [F b ( ®y * V + F c ] A dS (l8) 

while the pressure force is obtained from Eq. ( 17 ) to be. 



(19) 


This pressure force is combined with the shear stress integrated over 
surface ( 3 ) to obtain the total drag force (D^) on the jet flow due to 
blockage of the free stream flow. It is noted that for many injection 
situations the centrifugal force is in the opposite direction from 
the drag force, the drag force being in the positive n direction. 

Because of the complexity of the interaction between the jet and 
free-stream flows, the force resulting from the blockage effect is 
postulated to be the drag on an equivalent "solid" cylindrical jet 
shape inclined at an angle to the free-stream flow (after Abramovich 1 ^). 
This can be expressed as, 


D 


n 


IS 


C D ^ S 
n n 


n 


( 20 ) 


where the dynamic pressure of the free-stream flow perpendicular to 
the jet axis is given by 
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1 .,2 , . 

<3*0 = “2 P « V oo (21) 

n n 


Employing the definition V = V (e • I ) results in 

co 00 x n 
n 

- q » (e x e n 
n 


( 22 ) 


The reference area for calculating the drag force is the frontal area 
of the control volume (the surface facing the free— stream direction): 


s n * h. ds (23) 

where h is the local width of the jet measured in the t direction. 

Incorporating the expressions for q^ and S into Eq. (20) allows 

n n 

the drag force to he finally expressed as , 

D n “ C D ’ ®n )2 h ds (2*0 

n 

The third term in Eq. (10) accounts for the net influx of 
momentum into the control volume, which for the n direction, is due 
solely to the flux across the slanted surface of the cylinder. It 
is represented by the rate at which mass enters the sides of the 
control volume, see equation (3), multiplied by the free-stream vel- 
ocity component in the n-direction. Thus, 
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Jk*- dS) = \ h 

Equations (l8), (2k) > and (25) are combined to yield the conservation 
of momentum in the n-direction: 


II 




ds = 


V E ds 
00 
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(25) 
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(26) 


where the mass of the control volume is taken to be m = pAds. The 
expressions for the dot products in terms of the cartesian 
coordinates as the dependent variables and the distance along the 
trajectory as the independent variable are given in Appendix B, 
Substituting the appropriate definitions the n-momentum equation 
becomes , 


OAV 2 


®A( Poo - P) R H + C D ^ 

ds n 


& 


d 2 x 

+ E R S-| 

ds 


(27) 


where the radius of curvature (R) of the trajectory is defined in 
Appendix B. Previous studies investigating the 2-D injection 
problem (usually vertical injection) proceed to put this equation 
into a form with the slope of the trajectory 9 or the angular orien- 
tation, as the dependent variable. Since the vertical and lateral 
injection problems will be treated here as special cases of the more 
general 3-D injection situation, it is desirable not to specialize the 
equations in that way. The reader will see in Appendix C how this 
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equation is nondimensionalized and put into the form used in obtain- 
ing a numerical solution on the computer., The direction cosines of the 

jet velocity vector (i.e. e ) were chosen as the dependent variables 

s 

for several reasons , one of which is that the order of the governing 
eqns . is reduced by one. In addition, the algebra is thereby kept to 
a manageable level. 

It is appropriate to mention that at this stage of the develop- 
ment , the approach used in a number of previous studies"^ >^3-25 
been to assume the area growth of the jet along the trajectory, the 

rate of growth being based on data measurements where s/d. < 10. This 

1 — 

was done in lieu of solving the s— momentum equation. The area growth 
can be obtained by assigning a certain shape for the jet cross- 
section (e.g. , circular or elliptical) and by allowing the jet width 

pQ 

to grow at a. specified rate. Schetz and Billig used the expression 
for mass flow in the jet, that is, the continuity equation, to 
eliminate velocity in the n-momentum equation, the resulting 
expression then being integrated to obtain a solution for the jet 
trajectory. Typical results obtained by this procedure are presented 
in Fig. 4 for a jet with an elliptical cross-sectional shape and are 
compared with experimental data acquired from the photographs in 
Ref. 25- 

Comparison of the assumed cross-sectional areas with the values 
obtained from experiments shows that the two are in reasonable 
agreement in the proximity of the jet exit, but that the values 
diverge as the jet proceeds downstream, (It is noted that the 
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investigators who made the area-growth assumption were predominantly 
interested in the jet’s trajectory in the proximity of the injection 
point.) the measured areas indicating a much more rapid rate of jet 
growth than the assumed values . This trend is reflected in the tra- 
jectory information where good agreement between the predicted and 
experimental trajectories is noted in the initial region after jet 
injection, but poorer agreement occurs farther downstream. Another 
effect of assuming area growth is seen in the erroneous trends for the 
theoretical jet velocity, illustrated in the figure by the velocity 
deficit curve. As noted, the jet velocity begins to increase at 
some point along the trajectory. The reader is aware, of course, 
that as long as the jet injection velocity is greater than the free- 
stream velocity, the Jet velocity will decrease continuously along the 
trajectory and eventually approach the free-stream velocity value for 
downstream. 

From these remarks it is obvious that an alternative approach 
should be considered such that the Jet cross-sectional area is 
permitted to be an unknown in the governing equations . In order to 
do this it is necessary to have another equation to solve along with 
the continuity and n-momentum equations . The equation expressing 
conservation of momentum along the trajectory satisfies this need. By 
using this additional momentum equation, a more natural description 
of the Jet flow properties is obtained as illustrated in Figure 4. 

Although the area assumption is to be discarded, it is still 
necessary to provide information concerning the width of the jet in 
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order to calculate drag and shear stress terms in the governing 

19 

equations. One approach is to follow Abramovich and assume the 
growth of the jet width along the trajectory using an empirical 
expression based on limited experimental data. Another approach is 
to specify a shape for the jet cross section, and use this with the 
computed area to calculate the Jet width. This latter approach is 
more appealing because it is easier to Justify its use on the basis of 

16 

available experimental data. Keffer and Baines x ,for example, have 
shown that a jet initially having a circular cross-section transforms 
to a "kidney" shape as the Jet penetrates into the cross flow. This 
shape remains approximately the same with increase in s as 
illustrated in Fig. 5. 

Prior attempts by researchers at approximating the jet cross- 

section have been limited to the elliptical shape, where the circle is 

29 

a degenerate case. Hirsh , for example, approximated the jet cross 
section as a circle, which was a useful assumption in his development 
because then he assumed the Jet flow to be axisymmetric. It would 
appear that the ellipse would be a more suitable approximation, 
particularly near the injection point where the jet flow is deformed 
by the large pressure and shear stress fields. An ellipse with a 
major-to-minor axis ratio of 5?1 was employed in Ref. 19, and later in 
Refs. 23-25 , while Ref. 26 assumed a ^:1 ellipse. The approximation 
used in Ref. 26 had the added advantage of accounting for the change 
from the circular shape at the jet exit to the h:l ellipse at a 
specified point along the trajectory. 



32 



Figure 5.- Cross-sectional pressure contours of a jet injected nor- 
mally into a cross flow, Vi/Y^ = 2.2 (after Abramovich^ ) ; 
solid and dashed lines represent lines of constant total 
and static pressure, respectively, and the shaded areas 
indicate the potential core region. 
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For ilit* pretent, study , the Jet croap-noet ton nh n.po In rrruttick'I to 
be elliptical with a ratio of major-to-minor axes of !> to 1 , the major- 
axis being the Jet width. It is also assumed that at the injection 
point (s = 0) the elliptical area is equal to that of a circle with 
diameter d^. The change of Jet width, h, along the trajectory is 
accounted for by the expression. 



The change in h with distance along the trajectory was calculated 
with the present theory and is shown in Fig. 6 compared to experi- 

1 O 

mental data obtained from Kamotani and Greber . The data points 

from Ref. 18 were obtained from contours of constant velocity normal to 

a Jet cross section, where the velocity excess had decreased to 10$ 

of the maximum excess. The theory predicts the trend of increase in 

h/d i with increase in s/d i# although the predicted values are higher 

than the measured values at large s/d^. The small effect of VR on 

h/d^ near the injection point is reflected by the theory. The 

19 

empirical expression used by Abramovich in his theory is shown for 
comparison. 

The value of associated with the 5:1 elliptical shape is 

n 

taken to be 1.6 in keeping with the -equivalent "solid" body argument 

and is assumed to be independent of the Reynolds number of the flow 

26 

over the ellipse (Re^). Wooler et. al. used a value of 1.8 in their 

19 

analysis, while Abramovich used 3.0, a value which was pointed out by 
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23 

Schetz and Billig as "being totally unrealistic. An indication of 

how sensitive the theory is to the choice of cross-section shape and 

blockage coefficient can be seen in Figure 7 where experimental and 

theoretical trajectory results are compared for VR = 1*.0. The effect 

of changing ellipse-axis ratio is presented in part (a) where =1.6 

n 

is used in the theoretical calculations, while the effect of varying 

Cp Tor a specific ellipse-axis ratio (5:1) is shown in part (b). 
n 


As noted, increasing either ellipse-axis ratio or blockage coefficient 

results in progressively lower theoretical trajectories. The trends 

discussed here for the case of VR = k.O are typical of other injection 

velocities. The fact that the theoretical results obtained with 

C D = 1.8 (Fig. 7(b)) are in better agreement with the experimental 
n 

trajectories than the results obtained with C D = 1.6 is misleading. 

n 

It will be demonstrated later (Fig. kj) that the theory estimates 

mass flows in the Jet that are too low, so that if this disparity 

was corrected the results obtained with =1.6 would show improved 

n 


agreement with the experiment. 


s -momentum 

The s-momentum equation is obtained by taking the s-components of 
the various vector quantities in Eq. (10). The resulting expression 
represents a balance between the rate of change of Jet momentum and 
the forces on the jet due to changes in mean jet pressure, to buoyancy 




Figure 7«- Effect of Jet cross-sectional shape and blockage coefficient 
on theoretical trajectories for VK = k-, a = 90° 8 =0° 
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(b) 5:1 ellipse. 


Figure 7.- Concluded. 
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caused by density differences between the Jet and free-stream fluids, 
to entrainment of ambient fluid into the jet, and to shear stress 
between the jet and free-stream fluids. 

In the s— direction the shear and pressure forces are not combined 
as they were in the n— direction. The appropriate pressure force 
acquired from Eq. ( 17 ) is. 



where the local pressures are integrated over the control surface to 
get. 


F * - R 
Pc 


© A © + P © A ©~ P @ A ® 


(30) 


A Taylor expansion can be performed to obtain the pressure and area at 
surface (2) in terms of the appropriate variables at surface (T) , while 
PQ is taken to be the average of pq and , and Aq to be the 

difference between A^^ and Substituting these expressions into 

Eq. (30) and neglecting terms having higher orders of ds, the pressure 
force in the s-direction is found to be. 


F = A ds 
p s 3s 


(31) 


The force contribution of the surface stress tensor in the s-direction 
is approximated by* 
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h dS 


T 



d6 = TrhT ds 


(32) 


where t ^ is the local shearing stress in the jet flow acting on the 
slanted surface of the control volume. The integrated local shear 
stress is represented by an average shear stress t, which is assumed 
to act on a surface having the area irh ds. Although the physical 
problem is much more complicated than is described here, this approach 
will provide for a viscous force in the s-momentum equation, which has 
not been done in most of the previous studies using integral techniques. 

The contribution of the net influx of momentum into the control 
volume to the force balance in the s-direction is obtained by performing 
the operations suggested by the third term in equation (10) and is 
found to be : 
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( 33 ) 


The definition of averaged Jet properties was used to derive this 
equation and the momentum flux entering the sloping surface of the 
control volume was represented by the rate that mass flows across the 
surface multiplied by the free-stream velocity component in the 
s-direction. Equations (31), (32), and (33) are combined with the body 
force term from Eq. (l6) and the dot product expressions from Appendix B 
to obtain the s-momentum equation. 



ho 


= SA(P„ - P) g - A || - „hT ^ E g- (3M 

To evaluate the static pressure gradient along the trajectory, 

( 9p/3s) , the assumption is made that the free-stream static pressure 
field around the Jet perimeter imposes itself on the jet flow. This 
is the usual type of assumption made concerning other free turbulent 
processes, such as jet injection parallel to a mainstream (coaxial flow) 
or jet injection into a reservoir (free jet flow). For the present case 
where the jet structure is considered as an elliptical cylinder inclined 
at an angle to the free-stream flow, there are large variations in the 
free-stream pressure field around the Jet due to the blockage effect 
that the jet has on the free-stream flow. Some idea of the static 
pressure variation around the perimeter of a jet cross section, idea- 
lized as a circular cylinder, can be obtained by observing the experi- 
mental pressures in Fig. 8. An assumed pressure distribution to be 
used in the theory is also presented. As noted, the assumed pressures 
on the front of the cylinder (0 £ 6 £ n/2) are in functional form and 
were obtained from potential flow theory, while the pressures on the 
back of the cylinder (tr/2 < 0 £ it) are assumed to be equal to the 
free-stream pressure. Several researchers (e.g. Ramsey 17 and Kamotani 
and Greber ) have approximated the pressure field around the jet 
by examining the potential flow over various cylindrical shapes. 

Although the pressure field resulting from the turbulent jet injection 
process is very complicated and does not lend itself to be categorized 




1*2 


in this simple a fashion, the pressure variation in Fig. 8 is adequate 

for use in the present mathematical model. The local surface pressure 

(Co + p^) is used in the expression 
p n 

rTT 


P = 


1 

£ 


p d0 


(35) 


d0 


to obtain the average static pressure acting on the cylinder. Per- 
forming the integrations in Eq. (35) using the pressure distribution 
shown in Fig. 8 we obtain. 
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where it is recalled that q^ is the free-stream dynamic pressure 

n 

normal to the trajectory (see Eq. (22)). This equation implies that 

the average static pressure on the Jet cross section is less than the 

free-stream static pressure but approaches p^ as q^ approaches 

n 

zero. This occurs when approaches zero and/or when the Jet 

becomes parallel to the free-stream flow (i.e., a = 0). Assuming that 
the average pressure imposes itself on the jet flow (i.e., the 
pressure in the jet flow becomes p), p can be differentiated with 
respect to s to get: 
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It is noted that if the Jet trajectory iB in either the vertical (x-y) 
or lateral (x-z) planes, the above expression for the pressure gradient 
simplifies to, 


ff = (sin aHcos a) ~ (38) 
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which is the form used by Campbell and Schetz . 

One of the basic properties of a viscous fluid is that a velocity 
gradient or discontinuity in the flow results in a frictional shear 
stress. For the present case of a Jet injecting into a cross flow, 
we envision the shear stress in the s-direction to be proportional 
to the differences between the Jet velocity and the free-stream velocity 
component tangent to the Jet flow. These shear stress can be 
approximated by. 


T = p(v + e ) |~ (39) 

3U 

where U represents a velocity in the s-direction and — represents 

the gradient of that velocity in the n-direction. The kinematic 

viscosity (v) will be neglected for the present study since, for 

turbulent mixing flows, it is small compared to the virtual (eddy) 

viscosity, e. The method used to estimate the eddy viscosity for this 
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Jet injection process is based on Prandtl's hypothesis and is valid 
only for free turbulent flows. The viscosity is represented by. 
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where is an empirical constant and. b is the width of the mixing 

zone. The minimum velocity is defined as the free-stream velocity 

component tangent to the direction of the jet flow (V m ), while the 

s 

maximum velocity is defined as the mean jet velocity in order to be 
compatible with previous mean flow assumptions. This representation 
is illustrated in Fig. 9. 

The velocity gradient is approximated by. 


3U _ max min 
3n b 


(V - V ~) 
« ds 
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(41) 


so that the shear stress can be written as , 

T - ‘V 7 - v . fr> 


(42) 


Incorporating the pressure gradient term, Eq. (37), and the shear 
stress term, Eq. (42), into equation (34) yields the final form of 
the s-momentum equation: 
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Figure 9.- Representation of velocities used to estimate shear stress 
in s-direction. 


In order to obtain an estimate for the empirical constant K , it 
is necessary to rely on existing information related to less complex 
turbulent flows than the transverse jet injection process considered in 
the present study. The procedure used here is to estimate at the 

beginning of the injection by assuming that a two-dimensional free jet 
boundary exists between the jet velocity and the free-stream velocity 
component tangent to the jet flow. Schlichting's description of a 
free- jet boundary is used to obtain the expression, 

1 ^ «■*> 

where c is the rate of spread of the mixing zone and A is a velocity 
parameter. Although c might be expected to be a function of both s 
and V. /V , it is assumed to be a constant for the present study. The 
value of c is taken to be 0.32 which is more representative of the 
rate of spread observed for the experimental data in Fig. 6 than is 
Abramovich's value of 0.22. 

The velocities U and U . in the analysis of the free-Jet 
max min 0 

boundary are considered to be constant as the jet proceeds away from 
the point of initial flow interaction. In the present situation, 

U and U . continuously change as the jet flow is decelerating 
and bending over. Accordingly, the velocity parameter used in Ref. 

32 is re-defined as, 
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which forces to be dependent on the local velocity conditions along 

the trajectory. It is noted that at the initial point of a normal 
injection A = 1.0 which is the value Schlichting uses in his discuss 
sion. Equation (44) is used until attains the following value 

prescribed by a circular free-jet analysis : 



0.00217 

c 


(46) 


This description of is assumed to apply for the remainder of the 

jet's trajectory. An example of the variation of along the 

trajectory is shown in Figure 10 for several injection conditions. 

t-momentum 

The t-momentum equation is obtained by taking the t-components of 
the various vector quantities in Eq. (10). The resulting expression 
represents a balance of forces on the control volume due to buoyancy, 
to blockage of the free-stream flow, and to entrainment of ambient 
fluid into the jet. Similar to the derivation of the n-momentum 
equation, the pressure force is combined with the shear stress 
integrated over surface © to obtain the total drag force (D^) acting 
on the jet in the t-direction. Accordingly, 




Figure 10.- Variation of empirical constant Kj along the trajectory 
for several injection conditions; ou = 90° , g = 0°, 


where 


- «.<% • V (w: 

If the jet cross-section were circular, then would equal , 

"fc ft 

and S would equal S . However, since the shape is elliptic and 
"b ft 

not circular, then Ch ^ and S ^ S . Recalling that the 

t n h 

elliptical- shape has an axis ratio of 5 to 1, we have ds so 

that 
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where is taken to he 1.0. 

u 

The net flux of momentum in the t-direction entering through the 
sloping surface of the control volume is represented by the rate at 
which mass flows across the surface multiplied by the free-stream 
velocity component in the t— direction. This results in. 
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The body force term from Eq. (l 6) is combined with Eqs. (49) and (50) 
to obtain the t-momentum equation: 
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(51) 


After substituting for the dot products (Appendix B), the torsion (t ) 

associated with the trajectory is arranged so that it is in the 

numerator of the terms . The reason for this section is that x is 

o 

expected to have a small value for the present study (t - 0 for a 
two-dimensional trajectory). The resulting expression is. 
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For the case where the jet follows a two-dimensional path, the 
t-momentum equation is an identity (see Appendix C), and hence its use 
is not necessary in the procedure for obtaining a solution of the jet 
trajectory and flow properties. 

Heat- Energy 

Until this point only the mass and momentum aspects of the jet 
injection process have been discussed, however, since the present 
investigation is concerned with heated discharges, it is necessary 
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to also consider appropriate methods of describing the thermal 
characteristics of the flow. In particular, it is advantageous to 
determine the change in mean jet temperature resulting from the 
penetration of the jet into the cross flow. This can be accomplished 
by monitoring the heat loss from the control volume, this heat loss 
resulting from several heat-transfer mechanisms. 

The first type of heat-transfer mechanism pertains to the 
reduction in energy content per unit volume (pc p T) of the jet fluid due 
to the entrainment of free-stream fluid at a different energy level 
(pc T) . Applying this concept to the control volume results in the 

p 

expression , 


(m C p T) 2 = (m c p T) x + m e (c p T)^ (53) 

where (mc p T)^ represents the energy level in the control volume that 

would exist if there were no entrainment, and (me T) represents the 

P — 

equilibrium energy level resulting from the complete mixing of the jet 
and entrained fluids. Since the specific heat (c p ) of water is fairly 
insensitive to temperature changes, the various specific heats in 
Eq. (53) are assumed to have the same value. 

Forced convection, the second type of heat transfer mechanism 
being considered, results when the free stream flows around the heated 
jet fluid and extracts heat energy from the jet in the process. This 
heat transfer is analogous to the forced convection in separated flow 
over a heated cylinder, where the cylinder is cooled by the fluid 
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flowing normal to the cylinder's axis. To be consistent with our 
previous arguments, the convective heat transfer is estimated by 
considering the Jet structure as a cylinder inclined at an angle to 
the fre e-stream flow. 

33 

Eckert and Drake give several examples of film heat-transfer 
coefficients occurring in this type of flow situation and suggest the 
following expression for estimating an average Nusselt number: 

Nu d = 0.43 + 1.11 ^ ( Re d )^ (Pr) 0,31 (54) 

Values of Prandtl number for water at different temperatures were 
obtained from tables in Ref. 33 for use in the above expression. The 
Reynolds number is defined using the "effective" diameter of the Jet 
as the reference length and the free-stream velocity component 
perpendicular to the jet axis as the reference velocity, so that 
Nusselt number will be sensitive to the changes in local flow condi- 
tions as the Jet penetrates into the crossflow. The Reynolds number 
is thus , 
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For vertical or lateral injection this reduces to, 
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The values of Re^ occurring in the present study suggest the selection 
of ip = 0.45 and £ = 0.50 for use in Eq. (54). 

The definition of Nusselt number (hd/k) is used to obtain the 
average film heat-transfer coefficient (h) which yields the rate of 
heat loss from the jet fluid, 

Q = h BCT^ - T) (57) 

where B denotes the cylindrical area of the jet control volume. 

This, in turn, results in a temperature change in the jet flow due 
to this convective heat loss. 

An example of the temperature results obtained when these two 

heat-transfer mechanisms are incorporated into the analytical model 

is shown in Fig. 11, where the theoretical calculations were made 

with the same injection conditions as the data with VR = 5.2. The 

trend of temperature decrease along the trajectory measured in the 

present study is adequately estimated by the theory, the predicted 

average temperature values falling below the measured maxi m um 

temperatures . These results are substantiated by the temperature 

IT 

data for air injection processes measured by Ramsey and Kamotani 

1 8 

and Greber . The theoretical temperatures obtained by considering 
only the effects of entrainment are presented to demonstrate the 
relative magnitudes of the two types of heat-transfer mechanisms. 
Convection is seen to be the dominant mechanism for determining 
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Figure 11.- Variation of jet temperature with distance along trajectory 
for heated injection process; > T^. 
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temperature loss in the early stages of the jet injection process, 
while the effects of entrainment "become dominant as the jet proceeds 
downstream. 


Solution Procedure 

An iterative method is employed to obtain a numerical solution of 

the highly nonlinear governing differential equations at specific 

locations along the jet trajectory. Appendix C shows how the s,n, and 

t-momentum equations are non-dimensionalized and put into the forms 

used in the numerical technique. It was found that the equations could 

be simplified somewhat by using direction cosines (u and w) as the 

dependent variables rather than x,y, and z. This means that at each 

point (j+l) on the trajectory a solution to this initial value problem 

involves determining values for U J+1 » w j +1 » P_j +1 » A j+i» 811(1 V j+i' 

The basic solution procedure is to solve the s-momentum equation 

for the jet momentum in the control volume, where the coefficients in 

that equation are estimated using the flow property values obtained 

from the solution at the previous location on the trajectory. The 

s-momentum is used in conjunction with the continuity equation to 

provide. an update on V^ +1 and the heat loss from the control volume 

is calculated to provide new estimates for T.,, and o. . The 

j+l j+l 

current flow property values are then used in the coefficients of the 
n-momentum equation to obtain a solution for (du/ds) , where a 
central difference scheme provides Uj + g. Information acquired from 
the s- and n-momentum equations is used to solve the t-momentum 
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equation for (d w/ds ) , from which w 

J <J 

central difference scheme. The most recently calculated values of 
the direction cosines and flow properties are used to iterate back 
through the governing equations, convergence to a satisfactory solu- 
tion occurring in only a few iterations. 

Incremental values of x, y, and z are obtained from the final 
value of trajectory slope and from the assigned value for As. These 
increments are added to the coordinates of the previous location on 
the trajectory to obtain new x, y, z trajectory coordinates. This 
procedure is repeated at each incremental "step" along the trajectory 
to provide a solution for the trajectory and cross-sectional area of 
the jet, as well as the jet flow properties of mass, velocity, 
momentum, and temperature. The theoretical results presented herein 
were obtained with a constant incremental step size of 0.01 d^. 

For the case where the jet path is two-dimensional the 
information from the t-momentum equation is redundant; therefore, it 
is only necessary to solve the s- and n-momentum equations, along 
with Continuity, to obtain the desired solution. This is amplified in 
Appendix C where further observations are made concerning the 
numerical solution of the governing equations . 



IX . EXPERIMENT 


The experimental portion of the present investigation was con- 
ducted in a water channel located in the Hydraulics Laboratory at 
Virginia Polytechnic Institute and State University. A description of 
the test hardware and the test results is presented in the following 
section. 


Apparatus 

Water supply 

A centrifugal pump (1200 gal. /min. capacity) was the primary 
source for providing the desired flow conditions in the water 
channel. The sketch in Figure 12 shows the path of the water after it 
has been extracted from the sump by the pump, the bulk of the water 
being directed through the 10" line back to the sump. The necessary 
water flow for the channel was bleed from the 10" line by a 2" supply 
line. A globe valve located on the 2" line upstream from the water 
channel was used to control the flow rate to the channel. The actual 
flow rate before and after a test sequence of the main stream was 
determined with a weigh tank and a stopwatch. 

After the pump was started, the sump valve was adjusted to give a 
predetermined reading on a flow meter located on the 10" line upstream 
from the sump valve. This reading was used for all the tests using 
the pump so that repeating test conditions in the water channel would 
be possible. The centrifugal pump experienced problems with a foot 
valve after Run #37 (see test log in Chapter XIII ), so that water flow 
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in the channel was provided by city water for the remainder of the 
tests. The water from both sources was at room temperature and was 
clear and bubble free, as the photographs to be presented subsequently 
will show. A flow diverter was provided at the discharge from the 
channel so that water could be routed either into the weigh tank and 
hence into the sump, or into a drain. This diverter facilitated the 
flow rate measurements and provided a way to prevent the returning of 
dyed water (resulting from the jet injection process) to the sump. 

Main-flow channel 

The main- flow channel was a Hydraulic Demonstration Channel, 

Serial 116221, by Hydraulic Design and Products Company, Minneapolis, 

/ 

Minnesota and was constructed of plexiglass throughout except for some 
of the support points which were made of aluminum. The passage of the 
channel consisted of a 6" x 12" cross section approximately 8 feet 
long (Figure 13). Because of influences on the flow near the entrance 
to and exit from the passage, the usable test area was restricted to 
the region midway between the entrance and exit stations. The ch ann el 
was equipped with a sluice gate and weir both of which could be 
adjusted to alter the flow through the passage. The sluice gate was 
set so that the bottom of the gate was below the water level in the 
header, thus minimizing the number of air bubbles that passed through 
the gate to the channel's passage. The weir was adjusted to get the 
desired water depth at the jet exit stations. 
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Several arrangements were investigated for the discharging of the 
water from the 2" supply line into the header. Initially, a straight 
section of pipe was used to discharge the water into the header, but 
this arrangement resulted in very large disturbances in the flow in 
the test area of the passage. It was found that utilizing an elbow 
on the supply line so that the water is discharged towards the back of 
the header greatly reduced the disturbances noted in the channel. 

This is the arrangement shown in Figure 13 and was used throughout the 
tests . 

The basic water channel was altered in several ways for use in 
the present experimental study. Flow visualization tests indicated 
that , even with the elbow on the supply line , some degree of flow 
angularity and disturbances existed in the vicinity of the jet exit 
station. A flow straightening system was installed just downstream of 
the sluice gate to alleviate this problem. The system consisted of 
several layers of 1/8 inch wire mesh followed by a unit which held a 
bank of plastic soda straws { l/k inch in diameter). These straws were 
3 inches long, arranged streamwise, and were retained fore and aft by 
window screening. Subsequent flow visualization tests indicated that 
this flow straightening system provided a flow near the jet exit 
station that had few disturbances and had no apparent flow singularity. 

The second alteration was to drill 2.50-inch diameter holes through 
the bottom and one side of the channel 27 inches downstream of the flow 
straightener to permit insertion of the injection chambers. The third 
alteration was to cover the bottom and one side of the channel with 1/8 
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inch thick white plastic sheets that vere ruled with 1-inch squares 
in black paint. Each of the lines in these squares was 1/16 inch in 
width. 

Last , the channel was fitted with two instream copper-const ant an 
thermocouples located upstream of the flow straightener but downstream 
of the sluice gate. One of these was used to measure the main flow, or 
free-stream, temperature while the second was used in conjunction with 
thermocouples in the injection chambers to permit direct measurement of 
the temperature difference between the injected and free stream fluids 
(AT.). 

Injection system 

The various components of the injection system which control the 
flow rate and temperature of the injected fluid are shown in Figure lU . 
The injectant was water which was dyed with equal parts of red and 
green food coloring, 2 ounces per gallon. Propylene Glycol which has a 
specific weight of 6U - 8 lbf/ft 3 at 68°F, accounted for of the food 
coloring volume. 

An aluminum heating vessel 8" in diameter and 20" in length served 
as a reservoir for the injectant and had a fluid capacity slightly 
greater than 2-1/2 gallons. The fluid was driven from the heating 
vessel by means of air pressure (10 - 15 psig) which was controlled by 
a regulator. The water- flow rate was controlled by a needle valve and 
monitored with a Fisher-Porter P4 105711 flow meter. The tubing 
connecting the heating vessel to the flow meter and hence to the 




Figure 14.- Schematic of injection system. 
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injection chamber was l/k inch (o.d.) polyethylene to reduce heat loss. 
Heating of the injectant was accomplished electrically by a General. 
Electric 150-watt immersion heater that was regulated by a Variac. 

Fiber glass insulation material was used to reduce heat loss from the 
vessel. 

The injection chambers , machined out of solid lucite to minimize 
heat losses , were inserted through the channel wall and floor such that 
the surface facing the flow was flush with the white plastic surfaces. 
This is illustrated in Figure 15 which shows details of the wall 
injection chamber. A 9/1 6 " flat bottomed hole was drilled to within 
1/8" of the surface facing the flow, and a .0625" hole bored through to 
form the injection port. Two injection ports were used in the wall 
chamber during the duration of the experimental study; an injection port 
parallel to the chamber center-line was used for the lateral injection 
tests , and a port at a 45° angle to the center-line was used for the 
oblique injection tests. Both of these ports were located away from 
the chamber center-line so that a change in injection conditions could 
be obtained simply by rotating the injection chamber. For example, a 
change in vertical position for both ports was obtained by rotating the 
chamber, while an additional change in injection orientation occurred 
for the 45° port. The details of the injection chamber located in the 
channel floor are similar to those described for the wall chamber, 
except that the floor chamber had one 0.0595" injection port which was 
located on and parallel to the chamber center-line. 
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Figure 15.- Details of wall injection chamber used for lateral injection 
tests . 
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The outer portion of the 9/l6" hole was tapped for lA" NPT for 
connection with the feed tube by means of a compression fitting. A 
1/8" stainless-steel tube housing a copper-constantan thermocouple was 
inserted at approximately 1+5° from the outside into the 9/l6 inch hole. 
This thermocouple measured the injectant temperature and was positioned 
so that it was about 1 A" from the injection port. It was connected 
with one of the thermocouples in the free-stream flow to directly read 
the temperature difference between the injectant and the free-stream 
fluid. All of the thermocouples were coated with a thin epoxy coating 
to avoid corrosion, and although the epoxy coating increased the 
response times, calibration showed the response to be more than 
sufficient for the type of tests conducted herein. 

Instrumentation 

The two thermocouples placed in the free-stream flow were encased 
in 1/8" stainless-steel tubes , similar to the arrangement used to 
measure injectant temperature,. The outputs from these thermocouples 
were recorded on a Hewlett-Packard 7100B strip chart recorder. An 
electric ice junction by Joseph Kaye & Company was used to provide a 
reference for one of the thermocouples measuring free-stream temperature. 
The other thermocouple in the free-stream flow was connected with the 
thermocouple in the injection chamber by connecting the const antan 
leads of the two thermocouples together, and the respective copper leads 
to the recorder terminals. This arrangement allowed the direct 
measurement of the difference between the injected fluid and the free- 
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stream fluid. The millivolt output was converted to °F using a standard 
copper-constantan thermocouple calibration. Before each test sequence 
the outputs from the thermocouples were checked by using a beaker of 
water and an accurate thermometer. 

Two Graflex 4— in. x 5 - in. cameras were arranged above and to the 
side of the channel opposite the injection ports (see Fig. l6a) and 
were fitted with 127 mm, f/4.7, lenses. These lenses were operated at 
f/22, 1/400 sec to provide most of the photographic information obtained 
during the tests . Several photographs were obtained using a longer 
exposure (f/32, 1/15 sec) to provide a time-average of the fluid 
motion. For the free-stream velocities of this study, a fluid particle 
would move 3 to 6 jet diameters downstream during the long exposure. 
Illumination was accomplished with front lighting from a single 500-watt 
photo hood. Poloroid type 57, ASA 3000, film in sheet film holders was 
used by both cameras. 

The requirement that the pictures be taken simultaneously by both 
cameras lead to the setup illustrated in Figure 16(b). Heavy duty 
6-volt lantern batteries were used to provide the necessary current to 
a solenoid mounted on the side camera. By tripping the shutter on the 
top camera, a circuit is completed allowing the solenoid to trip the 
shutter on the side camera. 

A calibration of the Fisher-Porter flow meter was performed because 
the flow meter was marked to read pounds per minute for a fluid with ' 
specific gravity of 1.32, and the injectant for the present study was 
predominantly water having a specific gravity of 1.0. The calibration 
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(a) Camera arrangement. 
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<b) Shutter arrangement. 


Figure 1 6.- Schematic of camera arrangement and shutter operation. 
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was conducted for both the unheated and heated conditions of the 
injectant for a range of flow rates. The procedure was to adjust the 
needle valve to give a constant reading on the flow meter for a fixed 
period of time (usually 60 seconds), while collecting the injectant 
for a later temperature and weight measurement. The unheated injectant 
temperatures were about 70°F, and the temperatures of the heated 
injectant ranged from 125 F to l48 F. The results of the calibration 
are presented in Figure 17 where the measured values for flow rate of 
the heated injectant are represented by the ticked symbols. Since no 
discernible trend of temperature on flow rate was observed, a curve 
was faired through the data to represent the calibration for all 
temperatures. This calibration curve was used for all test conditions 
and provided the values of nu listed in the test conditions presented 
in Appendix A. 


Tests 

Experimental procedure 

The procedure for obtaining the temperature and photographic 
information for a given injection test began with the establishment of 
the desired flow conditions in the water channel. The centrifugal 
pump was used to provide the necessary flow in the water channel for 
the lateral injection tests, while city water was used for the channel 
flow for the vertical and oblique injection tests. Several measurements 
were made of the flow rate in the water channel using the weigh tank, 
and the weir was adjusted to set the water depth of the main flow at the 




Figure 17.- Calibration of flow meter for heated and unheated injectant. 
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injection station at a certain value. This water depth was 2-15/16" for 
the lateral injection tests, and 4—13/16" for the vertical and oblique 
injection tests . 

The injection system was checked to see that the heating vessel 
was supplied with an adequate amount of dyed injectant, and then the 
vessel was pressurized. If a heated jet test was to be conducted, 
approximately 15 minutes was necessary for the injectant to reach a 
desired temperature. The jet flow was turned on briefly to see that 
the flow meter and jet orifice were free from any restrictions, as well 
as to focus both cameras. When a test was ready to be conducted, the 
jet flow meter was set at a desired setting. If it was a heated injec- 
tant test, the needle valve was opened wide to ’purge’ the injection 
system with the heated fluid, before setting it at the desired setting. 
The effects of increased injectant temperature on the jet flow charac- 
teristics were investigated by setting the flow rate of the jet at a 
fixed value regardless of the temperature of the injected fluid. As 
is shown in the test log in Appendix A, the resulting injection veloci- 
ties are practically the same since the density of the injectant has only 
a small variation with temperature within the temperature range 
investigated. 

When the photographs were taken simultaneously with the two 

cameras the strip chart was marked so that the temperature readings 

would coincide with the photographs. The strip chart recorded continuous 

readings of the two temperature measurements. The values of T and 

00 

AT^ shown in the test log were those values read off the strip chart at 
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this point. 

Test conditions 

The tests were conducted for the submerged injection of the jet 
fluid into the main flow in three ways . These will be referred to as 
lateral, vertical, and oblique injection throughout the report. A 
cartesian coordinate system is defined such that the x axis is 
parallel to the direction of the free stream flow, and the y axis is 
in the vertical direction (i.e. in line with the gravity force). The 
z axis is the perpendicular to both the x and y axes. The orientation 
of the jet's axis at the injection point is defined by using the angles 
ol and 9 where a is the angle between the jet axis and the x 
axis, and 3 is the angle between the jet axis and the y axis. The 
angular orientation of the jet axis is shown in Table II for the 
different injection situations, the jet fluid being injected perpendicu- 
larly into the main stream for the lateral and vertical injection cases. 
The vertical distance of the injection port from the channel floor (a) 
and the horizontal distance from the flow straightener (L) are defined 
in Figure 18 to aid in describing the locations of the ports in the 
water channel. These distances are listed in the test conditions that 
are tabulated in Appendix A. 

The different orifice locations for the lateral and oblique 
injection cases were obtained by rotating the wall injection chamber. 

For lateral injection, this gave the effect of changing the distance 
from the jet flow to the free surface of the main flow. For oblique 
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injection, this completely changed the orientation of the injected 
jet with respect to the main flow, as evidenced by the ol and 3. 
values shown in Table II. Only one orifice location was investigated 
for the vertical injection case and it was located on the channel floor 
mid-way between the two channel walls . 


TABLE II.- Angle of jet axis at injection point. 


Injection 

a. 

i 


Orifice 

Lateral 

90° 

NO 

o 

o 

1 and 2 

Vertical 

o 

o 

ON 

0° 

3 

Oblique 

o 

o 

ON 

50° 

k 


H 

s=- 

o 

o 

90° 

5 


40° 

o 

o 

ON 

6 


The experimental tests were conducted for the injectant having a 
range of injection velocities and temperatures, and the resulting test 
conditions are presented in tabular form in Appendix A. Each test is 
assigned a Run number for easy identification with the experimental 
results presented in the text. The tables list the measured values 
for the mass flows of the jet and free-stream flows, and for the 
temperatures of the jet and free-stream fluids. An explanation of how 
the other properties are deduced from the measurements is included in 
the appendix. The reader should note that the values for velocity ratio 
(VR) and temperature difference (AT ) presented in the text are 
"rounded off" versions of the values listed in the tables. 
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(a) Lateral injection. 
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(b) Vertical injection. 
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(c) Oblique injection. 

Figure 18.- Jet orifice location for different injection conditions; 
a and L are shown in inches. 
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The physical dimensions of and flow properties in the water channel 
were such that fully developed flow (where the mean velocity profiles 
across the channel are not changing with distance downstream) could not 
occur. Using nominal flow conditions for the vertical and oblique 
injection tests ( = 0.33 ft/sec and = 68.0°F), the Reynolds 
number based on the length along the channel measured from the flow 
straightener was 6.8 x 10^ at the injection station. Similarly, Re^ 
was 12.8 x 10^ for the lateral injection tests. These values of 
Reynolds number imply a laminar condition for the free-stream flow. 
Assuming that the flow along the channel walls is equivalent to laminar 
flow along a flat plate, rough estimates for the laminar boundary-layer 
thicknesses on the channel walls at the injection point were found from 
equation (58) to be 


- 5 — 
L ssr 

•Li 


(58) 


0.5 inches for the vertical and oblique injection tests and 0.38 inches 

for the lateral injection tests. In terms of jet diameters, 6/d^ would 

be 8 and 6, respectively. These values for 6/d^ are much larger than 

those usually encountered for air injection tests. For example, the 

IL6 37 

studies of Kamotani and Greber and Ramsey had <5/d^ values of 0.1 
and 0.6, respectively. Because the estimated boundary-layer thicknesses 
are many times larger than the injection port diameters, it is natural 
to express concern about the influence of the nonuniform free-stream 
velocity field on the injected jet characteristics. Since no 
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experimental measurements were made of free-stream velocities, an 
exercise was performed to assist in a better definition of the free- 
stream environment into which the jets were injected. This information 
is presented in Appendix A. 

It should be mentioned that the laminar condition of the main flow 
and of some of the inject ant flows (to be shown subsequently) is not 
compatible with the turbulent nature of actual, full-scale problems. 
However, since the flow in the mixing region of the jet became turbulent 
at or near the injection point in all test cases, it is believed that 
the simulation is reasonable . Also , no attempt was made during the 
tests to predetermine the relative magnitudes of the buoyancy and 
momentum forces in the jet flow. The ratio of inertia to buoyancy 
forces at the injection point (Froude number, Fr^) is presented in 
Appendix A for the variety of tests conducted. 
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RESULTS 

The three-dimensional trajectory information for the different 
injection orientations was obtained for the jet flow having a range of 
initial velocities and temperatures. The primary output of these 
experimental tests is photographs which are presented in this section. 

In all of the photographs the free stream is flowing from right to 
left. The "effective" velocity ratio (VR) utilized thro ugh out the 
present paper has been suggested by Margason 13 and Kamotani and Greber 18 
as the proper parameter with which to compare different injection 
situations . 

Lateral Injection 

Tests were first conducted with the injectant and free-stream 
temperatures equal in order to provide a basis of comparison for the 
heated jet results. Photographs of the unheated lateral injection 
process are presented in Fig. 19 for three injection velocities and 
show that an increase in Jet injection velocity results in further 
penetration of the jet into the mainstream. The photographs also 
indicate that for the higher injection velocities the area occupied 
by the Jet fluid begins to grow immediately after injection , and this 
growth continues as the jet structure bends over under the influence 
of the free-stream flow and body forces. At some point downstream of 
the injection station (dependent on the magnitude of the jet velocity 
compared to the free-stream velocity), the jet flow becomes parallel 
to the free-stream flow. The turbulent nature of the injection process 
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Side view 



Top view 



(a) VR = 5.3, Run 9. 




Figure 19 •- Photographs of lateral 
injection velocities; 


injection process for a range of 

AT. = -2°F. 
i 
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can be seen by the very irregular boundary of the jet which is 

indicative of large-scale eddies in the flow. Abramovich 1 ^ and Keffer 
3 _ . 16 

and Baines provide excellent descriptions of the evolution of a jet 
injecting into a cross flow. 

The basic characteristics of the injection process are the same 

for the lowest injection velocity (VR = 5 . 3 ) except that the jet 

fluid appears to penetrate a discrete distance into the free-stream flow 

before beginning to spread. This indicates that the Jet flow is 

laminar at the injection point and undergoes transition to turbulent 

flow at some point along the trajectory. This laminar flow situation 

results because the Reynolds number of the jet flow at the injection 

point (Re ) is less than the critical value of Reynolds number 
i 

(-2300) below which the Jet flow cannot be turbulent. For example, 

Re^ * 1U00 for the VR * 5*3 case. The higher injection velocities in 
a i 

Fig. 19 result in larger values of Re and hence allow the jet flow 

a i 

to begin spreading immediately after injection. The injection studies 
of Kamotani and Greber and Hoult and Weil^ are concerned with jets 
whose flows are initially laminar. In Hoult and Weil's investigation a 
buoyant plume issuing from a smoke stack was experimentally simulated 
and several tests were performed to determine when the plume became 
turbulent . 

Centerline trajectories were measured from the photographs in 
Fig. 19 and are presented in Fig. 20 compared with trajectories 
estimated by the present theory. As noted, the theory gives a 
reasonable estimate of the path of the jet as it proceeds downstream 
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and predicts the further penetration of the jet into the mainstream 
that results from increased injection velocity. Unless otherwise noted, 
the average free-stream velocity (V^) listed in Appendix A is used in 
the theoretical calculations which are compared with the present 
experimental data. 

One of the interesting flow phenomena resulting from a fluid 
injecting into a cross flow is illustrated in Figure 21 which presents 
several long exposure photographs for lateral injection with VR = 5.3. 
The tendency for the Jet fluid to gather or collect away from the jet 
center line (side view) and towards the rear of the jet flow (top view) 
is indicative of the strength of the rotational velocity field, usually 
interpreted as a pair of counter-rotating vortices (Fig. 22). These 
vortices are a product of the complex interactions that take place 
between the jet and free-stream flows and have been discussed and 
measured by a number of researchers 1 ^™ 20 . Keffer and Baines 16 , in fact, 
suggest that far downstream the limiting condition for the jet injec- 
tion process is a pair of counter-rotating turbulent line vortices 
moving with the speed of the main stream. During the present study 
several observations were made with flow visualization devices to 
confirm the presence of the vortices, however no attempt was made to 
measure the vortex size or strength. 

The effect of adding heat to the injectant is illustrated by the 
photographs in Figure 23, where the injectant temperature is at least 
5^°F greater than the temperature of the free-stream fluid. The gross 
features of this heated jet injection process are the same as those 
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Run 13 


Figure 21.- Long exposure photographs of lateral injection process; 
VR = 5-3, AT^ = -1°F. 


c-a 
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Figure 22.- Diagram of the interaction resulting from jet injection 
into a cross flow (after Abramovich 1 ?) . 
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Side view 


Top view 




(b) VR = 17. 9, AT. = 65° F, Run 31. 


Figure 25.- Photographs of heated, lateral injection process for 
several injection velocities. 
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discussed for the unheated jet in Figure 19. One noticeable effect is 

that the fluid for the heated jet appears to spread at a slightly 

greater rate immediately after injection; this is particularly evident 

for the lowest injection velocity case (VR = 5.2). This effect on 

spreading is related to the larger values of Re 3 which result when 

d i 

the jet fluid is heated. 

Comparing the photograhs of Figures 19 and 23 we see that there 
is no discernible effect of the higher injectant temperatures on the 
jet trajectory. The lack of any significant effect of temperature on 
the trajectory is believed to be due to the combination of three 
factors: (l) water density is relatively insensitive to temperature 

change, (2) momentum forces dominate any buoyancy forces present at 
these conditions, i.e. the Froude number representing the flow is very 
large, and ( 3 ) the jet fluid experiences a rapid heat loss along the 
trajectory. To obtain an idea of how rapidly the jet flow losses heat 
jet temperature was measured with a thermocouple probe at several 
locations along the trajectory. Time averages of these temperatures 
are presented in Figure 11 for VR = 5.2 and are shown plotted in the 
form of a temperature deficit. The measurements indicate a rapid drop 
in jet temperature as the jet begins to penetrate into the free-stream 
flow, and that after this rapid drop the jet temperature slowly 
approaches the free-stream temperature (T^) with increased distance 
downstream. The photographs in Figure 2k were taken in conjunction 
with two of the temperature measurements for the VR = 5.2 case. 
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The photographs presented in Figure 25 illustrate the effect of 
decreasing the vertical distance between the injection port and the free 
surface of the water channel. The basic characteristics of the 
injection process for the three injection velocities shown are the same 
as those discussed in Figure 19 for the lower position of the injection 
port. The one obvious effect of injecting closer to the free surface 
is observed for the highest injection velocity (VR = 28.5). The wave 
pattern that is generated on the free surface (top view) is an 
indication of how severely the free-stream flow is disturbed, or 
blocked, by the jet flow. 

Vertical Injection 

The effect of increasing the velocity of a jet injecting vertically 
into the main stream is shown in the photographs of Figure 26 where the 
injectant and free-stream temperatures are essentially the same. The 
basic characteristics of the injection process are the same as those 
discussed previously for lateral injection. The jet flow for VR = 8.8 
shows evidence of laminar flow immediately after injection, similar 
to the lateral jet flow (VR = 5.3) noted in Figure 19. In fact, the 
initial laminar region appears to be more pronounced for the vertical 
injection case than for the lateral injection case. In comparing these 
two runs, it should be noted that the free-stream conditions are 
different. The vertical injection situation has a lower free-stream 
velocity than the lateral injection case, and hence has a larger 
boundary layer thickness at the injection point. This larger value of 



88 


Side view 


Top view 



(a) VR = 5.0, Run 16. 



(b) VR = 17.2, Run 15. 



(c) VR = 28. 5, Run i4. 

Figure 25.- Photographs of lateral injection near the free surface for 
a range of injection velocities; AT\ = 0°F. 







Figure 26.- Photographs of vertical injection process for a range of 
injection velocities; AT = -3°F. 
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6 coupled with a lover value of Re, can account for the more pro- 

d i 

nounced laminar flow condition for the vertical injection than for the 
lateral injection. 

The center-line trajectories for the three injection velocities 
presented in Figure 26 are shown in Figure 27 compared to the theoreti- 
cal trajectories (solid line) which represent fully turbulent jet 
flows. As noted, the experimental mixed- flow jets penetrate much 
further into the free stream than would be theoretically expected for 
fully turbulent jets. In order to estimate the trajectories for this 
mixed- flow situation, the present theory was adjusted to account for 
the initial laminar portion of the jet flow. This was accomplished 

by assuming that the jet begins its turbulent growth at a point (y /d ) 

o i 

specified in the photographs in Fig. 26. Since the location and extent 

of the transition region in the flow are functions of injection 

conditions (e.g. Re^ ), as well as free-stream conditions (e.g. VR) , 

i 

it is expected that the values of y Q /^ will change accordingly. 

The appropriate values of 7 Q /^ used to modify the theory are shown 
in Fig. 27 and the resulting calculations represented by the dashed 
lines . 

The effects of increasing injectant temperature on the vertical 
injection process are presented in the photographs in Figure 28. The 
slight increase in spreading of the jet fluid with increase in 
temperature is seen in the jet flow immediately after injection. This 
is particularly the case for the Jet with the lowest injection velocity, 
this trend being similar to the temperature effects discussed for 




Figure 27 „- Experimental and theoretical trajectories for the vertical 
injection process, AT. = -3°F; experimental trajectories 
are for jets with initially laminar flow. 
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Figure 27.= Continued. 
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Figure 27. ~ Concluded. 



9b 




(b) VR = 29. 0, ATj = 25°F, Run 46. 


Figure 28.- Effect of increasing injection temperature on vertical 
injection process for a range of injection velocities. 




Figure 28.— Concluded. 
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lateral injection. The photographs for vertical injection indicate 

that increasing injectant temperature tends to decrease the amo unt of 

penetration for the jet, even though heating the injectant causes 

slightly higher values of VR. This result is better illustrated in 

Figure 29 where the center-line trajectories are shown for the various 

injection velocities and temperatures investigated. The effect of 

injectant temperature on trajectory is mainly the result of the 

temperature effect on the transition of the jet flow, where an increase 

in temprature increases Re a and hence decreases the extent of the 

d. 

1 

laminar portion of the flow. Theoretical calculations are shown for 
comparison with the experimental trajectories obtained with the 
highest injectant temperature, where the fully turbulent flow condition 
is represented by a solid line, and the initially laminar flow condition 
by a dashed line. The procedure for adjusting the theory to account 
for the initial laminar flow is the same as that reported for the 
results in Fig. 27, although the values of y^/d^ are smaller. 

Oblique Injection 

The primary characteristics obtained by injecting obliquely 
(a. = 90°, B^ = 50°) into the main stream are seen in Figure 30 to be 
similar to those discussed for the lateral and vertical injection 
processes. The condition of laminar jet flow near the injection 
point that was observed for the lowest injection velocity cases for 
the previous injection processes is not as apparent here even though 
the value of Re^ is comparable. This may be due to the increased 
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Figure 29.- Effect of injectant temperature on trajectory for vertical 
injection process; experimental trajectories are for jets 
with initially laminar flow. 
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Figure 29.- Concluded, 



Figure 30.- Photographs of oblique injection process for a range of 
injection velocities j cl = 90°, = 50°. 
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vorticity interaction that occurs between the jet flow and the wall 
boundary layer. For this type of oblique injection the effects of 
increasing injection velocity on trajectory are observed in two 
planes instead of just one plane as is the case for lateral and 
vertical injections. Centerline trajectories in the X-Y and X-Z 
planes are measured from the photographs (Figure 30) and presented in 
Figure 31. From these data it is evident that the jet flow penetrates 
farther into the free stream in the lateral (Z) direction than in the 
vertical (Y) direction. This trend is more obvious at the higher 
injection velocities and results because the lateral momentum of the 
jet is larger than the vertical momentum at the injection point where 
= 50°. The projections of the trajectory onto the two planes 
would be expected to be identical if the oblique injection process 
had = 1+5° and was void of any wall or boundary-layer effects. 
Trajectories estimated by the present theory are shown in Fig. 31 
for comparison with the experimental data. The theory does predict 
the experimental trends discussed above, although it estimates slightly 
less vertical penetration and more lateral penetration for a given 
value of VR than is indicated by the data. 

The effect of adding heat to the inject ant can be seen in the 
photographs in Figure 32, where increasing injection temperature causes 
a slightly greater rate of spreading by the jet fluid near the 
injection point, particularly for the lower injection velocities. 

Also, the higher temperature jets generally penetrate farther into 
the free-stream flow than those in Figure 30. The effect of injection 
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Top view Side view 



(a) VR = 8. 1, ATj = 46° F , Run 60. 



(b) VR = 15.9, AT j = 44° F, Run 59. 

V 



(c) VR = 26. 7, AT j = 48°F, Run 58 


Figure 32.- Photographs of heated, oDlique injection process for a 
range of injection velocities; on = 90°, g = 50°. 



'R = 38.0 


, AT j = 63°: , Run b' 


Figure 32.- Concluded. 
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temperature on centerline trajectory for oblique injection is shown in 
Figure 33 for the different injection velocities investigated. 

The fact that the oblique injection condition demonstrates an 
increase in Jet flow penetration with increase in injectant 
temperature, while the vertical and lateral injection conditions do not 
demonstrate this effect, draws attention to a subtle difference between 
these injection cases which is related to the orientation of the counter- 
rotating vortices in the respective Jet flows. For the lateral and 
vertical injection situations, the jet is injected perpendicularly 
through the boundary layer on the adjacent wall. The vortices formed 
during these injection processes are located symmetrically on either 
side of the jet center-line, such that a line drawn between the vortex 
centers would be parallel to the vorticity vector associated with the 
boundary-layer velocity field. This description is not completely 
accurate for oblique injection through the boundary layer. From 
observations of the oblique injection experiments, it appeared as if 
the vortex pair was "twisted" immediately after injection so that the 
pair was oriented in a manner similar to the vortices generated by a 
lateral injection, that is, where the line between the vortex centers 
is vertical. Although detailed measurements are necessary to 
validate these comments, it is believed that the observed effects of 
injection temperature on the jet path for the three injection conditions 
are due to differences in vortex formation in the respective jet 


flows . 




s® 
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Figure 33 = =* Continued 
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Figure 33.- Concluded. 
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Photographs of the oblique injection in the upstream (ol = ll+0 o , 
$ i = 90°) and downstream (cl = Uo°, = 90°) directions are 
presented in Figs. 3*+ and 35, respectively. Trajectory data 
approximating the centerline of the Jet paths in the photographs are 
shown in Figs. 36 and 37* For the case where ol = 1^0° the Jet 
reaches a maximum distance upstream before it proceeds in the down- 
stream direction. The data of Platten and Keffer 11 * also illustrate 
this trend (Fig. 36), even though the injection conditions are 
somewhat different from the present experiment. The trajectory 

10 

estimated by Ivanov's empirical expression , Eq. (59), 



is presented for comparison with the present data. 

For the case where the Jet is injected downstream (Fig. 35) the 
Jet path has less curvature and penetration than was observed for 
lateral injection (Fig. 19). This is reasonable since the curvature 
of the trajectory is governed by the forces perpendicular to the Jet 
path, which are proportionately less for cl = Uo° than for a. = 90°. 
In the limit as 0^ -*■ 0°, these normal forces go to zero resulting 
in the familiar co-axial Jet flow process whose trajectory is the X- 
axis. The present theory adequately estimates the trajectories for 
injection with ol = 40° (Fig. 37). 



Run 67 


Figure 3^.- Photographs of oblique injection in upstream direction 
(a. = l40°, $. = 90°) for VR = 8.7 and AT. = 6T 0 F. 
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Side view 


Top view 



(b) VR = 28.5, ATi = 43°F. Run 71. 


Figure 35.- Photographs of oblique injection in downstream direction 
(cl = 1*0°, 3^ = 90°) for several injection velocities. 
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Figure 37-- Experimental 
ct. = 40° , 6. 


and theoretical trajectories 
= 90°. 


for oblique injection downstream; 


115 


X. FURTHER EXAMINATION OF THE THEORY 


The purpose of this section is to discuss some of the limitations 
of the theory developed in Chapter VIII, and to demonstrate its 
versatility for handling a variety of injection situations. In order to 
establish the authenticity of the present theoretical method for 
estimating jet flow properties , its predictions are compared to 
estimates from other analytical models as well as to experimental data 
acquired from a number of studies. The last portion of the section 
presents a theoretical example of a jet with a three-dimensional 
trajectory. 

Two-Dimensional Trajectory 

Experimental trajectory data obtained from different investigations 
of air jets are presented in Figures 38 and 39 and show the two- 
dimensional paths of the turbulent jet for a range of injection velo- 
cities and orientations. These data were obtained from hot-wire 
measurements and, thus, represent the path that is traced by the 
maximum velocity in the jet flow. The theoretical trajectories were 
calculated with the same injection conditions as the experimental data 
and are in good agreement with the measured trajectories throughout 
the range of injection velocities and orientations. 

Theoretical trajectories calculated with the present theory are 

compared with theoretical and experimented results of other researchers 

in Figures 40 - 1+3. The analytical methods of Abramovich 1 ^, Schetz and 
23 2l+ 

Billig , and Reilly provide jet trajectories which are comparable to 
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Figure 38.- Experimental and theoretical trajectories of an air jet 

having a range of injection velocities; a. = 90° „ 
ft = n° 






2 


4 


8 


10 


x/dj 

(c) VR = 8.32. 


Figure 39.- Concluded. 
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those of the present theory for the injection conditions presented 
(Fig. 40 and Hi). It is recalled that the theories of Refs. 19, 23, 
and 24 assumed the growth of the jet cross-sectional area along the 
trajectory using an empirical expression (see Fig. 6) based on 
experimental data for s/d^ <_ 10. As a consequence, the trajectories 
predicted by these theories agree quite well with experimental data in 
the vicinity of the injection point. Care must be exercised in using 
these theories to estimate jet trajectories and flow properties at 
large s/d^ values. 

One of the best theoretical methods prior to the present is 
29 

that of Hirst , who attempts to account for the complex flow 

processes that take place as the flow evolves from a momentum jet near 

the injection point to a buoyant plume at large distances downstream. 

His results are compared with the present theory in Figs. 42 and 43 

for a range of injection velocities and angles. As noted, the 

present theory is in better agreement with the bulk of experimental 

data for all of the injection conditions. Since Hirst assumed a 

Gaussian type of velocity distribution in the jet, his theory is 

applicable only in the region where the jet flow has become fully 

developed. This explains why his theoretical trajectories do not 

originate at the injection point. The experimental data obtained by 

Gordier are shown in Fig. 42 because Hirst compared his theory with 

these data in Ref. 29. Gordier 's data, however, indicate greater 

17 

penetration by the jet than is seen for the other data. Ramsey 
suggested that this discrepancy was probably due to injection into a 
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Figure 40.- Comparison of theoretical trajectories estimated by present 
theory and by Abramovich; <x = 90° 9 = 0°. 



x/dj 


Figure ^1„- Comparison of theoretical trajectories estimated by present 
theory with those of Schetz & Billig and Reilly ; a. = 9C? 9 



(a)VR* = 4 0 

Figure 42,- Comparison of theoretical trajectories estimated by present 
theory Q with those of Hirst for normal injection j a. ~ 90° 


Experiment 

Keffer and Baines 16 
Kamotani and Greber 1 ^ 
Jordinson 15 
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Figure h2.- Concluded. 



(a) VR = 4. 05. 


Figure 43.- 


Comparison of theoretical trajectories 


estimated 


present 


S> 
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cross flow with a very thick boundary layer. This trend will be shown 
later in this chapter. 

Jet Flow Properties 

Examples of some of the theoretical flow properties obtained in 
the process of solving the governing conservation equations are pre- 
sented in Figures 41+ to 49. It was observed in the experimental 
portion of the present study that the cross-sectional area of the jet 
continually increases as the Jet proceeds along the trajectory. This 
trend is shown in Fig. 44 where the Jet area, normalized by the Jet 

area at the injection point, is plotted as a function of s/d^ The 

1 ft 

data points acquired from Kamotani and Greber's work were obtained 
by measuring the area encompassed by a contour of Jet velocity where 
the velocity excess had decreased to 10$ of the maximum excess with 
respect to the free-stream velocity component tangent to the jet flow. 
It is noted that Keffer and Baines 1 ^ also defined the edge of the Jet 
flow in this fashion. Experimental areaB for the VR = 30 case were 
obtained from the photographic information in Fig. 19 by assuming the 
cross-sectional shape to be a 5:1 ellipse and by measuring the minor 
axis. The data in Fig. 44 indicate that an increase in injection 
velocity (VR) results in larger rates of area growth with s/d^. 

Theoretical areas are shown for comparison and predict the same 
trends as the experimental areas, however, the theory underestimates 
the magnitude of the area growth experienced by the jet for the range 
of injection velocities shown. It should be mentioned that these 



Figure 44,- Variation of jet cross-sectional area with distance along 
trajectory for a range of injection velocities, a = 90° 0 
(Values from Ref. 18 are measurements of the area^bounded 
by a contour of jet velocity). 
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theoretical estimates of jet area are very sensitive to the amount of 
entrainment, small increases in entrainment resulting in large 
increases in the cross-sectional areas. It will be shown shortly that 
the entrainment model used in the theory could be improved, which would 
result in better agreement between the theory and the experiment in 
Fig. 44. 

As the area occupied by the Jet fluid grows with increase in 

distance along the trajectory 9 the jet velocity correspondingly decays. 

This is illustrated in Fig. 45 where the jet velocity is nondimen- 

sionalized by the injection velocity and where the V /V values for 

00 jL 

the various injection conditions are depicted by the arrows. The 

experimental data that are shown are measurements of the maximum jet 

velocity for VR ® 4, 8, and 00 (free jet). These data indicate that 

the trend for velocity decay is similar for all injection velocities. 

For the VR = 4 case there is a short distance (potential core) 

where the maximum jet velocity remains equal to the injection velocity. 

Further increase in s/d i results in continued decreases in jet 

velocity which eventually approaches the free-straam velocity value 

(V /V ) . Increasing injection velocity increases the potential core 
00 1 

18 

length, a maximum value being obtained for the free jet s and 
decreases the value of V /V that the Jet velocity must approach 
(note arrows), V/V = 0 for the free jet. The combination of these 
experimental trends helps in comprehending the theoretical velocity 
decay curves presented in the figure for a range of injection velocities. 
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Figure 1*5.- Effect of injection velocity on jet velocity decay with 
increase in distance along trajectory; cl = 90 °. 
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Comparing the theoretical velocities with the experimental data for 
VR = 4 and 8, it is observed that the theory predicts V/V. values 
that are greater than the experimental, values and which approach 

at a slower rate. The fact that the mathematical model est ima tes 
average jet velocities explains why the theory predicts a decrease in 
velocity immediately after injection compared to the existence of a 
potential core demonstrated by the maximum velocity measurements. 

The effects of increased injection velocity on velocity decay, 
which were discussed in Fig. 45, are put into perspective in Fig. 46 
by plotting the velocity as a velocity deficit, where the difference 
between the jet and free-stream velocity iB divided by the difference 
at the injection point. Presenting the results in this fashion causes 
the velocity deficit to approach zero as V . As noted, an 

increase in VR results in a corresponding increase in experimental 
velocity deficit at a given distance on the trajectory. This trend 
also essentially applies to the theoretical velocity deficit variations. 
Keffer and Baines 1 ^ observed that velocity deficit showed a universal- 
ity (i.e. independent of VR) when plotted against the distance from 
the virtual source of the jet flow. 

The variation of theoretical mass flux in the jet with distance 
along the trajectory is shown in Fig. 47 for a range of injection 
velocities. As would be expected from our consideration of mass 
conservation in Chapter VIII, the mass flow in the jet increases with 
increased distance along the trajectory, the higher mass flows occurring 
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Figure 46=- Variation of jet velocity deficit with distance along 

trajectory for range of injection velocities; a = 90° 

i 
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Theory Experiment 

VR Ref 



Figure kj.- Variation of mass flow in jet with distance along trajec- 
tory for range of injection velocities; ol = 90°. 
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for the lower injection velocities. This increase in mass flow with 

increase in s/d ^ is demonstrated by experimental data 1 ® for VR = k , 

although the measured values indicate higher mass flows in the jet than 

are predicted by the theory. A further indication that the theoretical 

mass flows should be greater than they are is provided by the measure- 

22 

ments of mass flow in a free jet. The fact that these mass flows 

are less than the experimental values for the jet with VR = k supports 

the premise that the entrainment rate for a jet in a cross flow should 

be greater than that for a free Jet. At certain values of s/d the 

1 

disturbing situation exists that the theory predicts mass flows that 
are less than the free-Jet values. Improvements in the entrainment 
model used in the current analytical effort will provide improved 
estimates of the Jet's mass flow, and consequently will yield more 
realistic cross-sectional area and velocity decay results than were 
observed in the last few figures. In the study by Schetz and Billig 23 , 
the mass flow was assumed to remain constant at the initial value. 

Their assumption is represented by the horizontal line in the figure 
and becomes more realistic nearer to the injection point. 

The variation of theoretical jet momentum flux with distance 
along the trajectory is shown in Fig. U8 for a range of injection 
velocities. The jet momentum, normalized by the value at the injection 
point, is seen to decrease immediately after injection, reaching a 
point on the trajectory where it attains a minimum value, after which 
it increases for the remainder of the trajectory. This trend occurs 
because the component of momentum in the y- and x-directions decreases 
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Figure 1+8.- Variation of theoretical jet momentum with distance along 
trajectory for range of injection velocities j a. = 90° 

S. = 0 °. 1 



139 


and increases, respectively, along the trajectory of a jet that is 
injected normally into a free-stream flow. For this situation (a. =90°) 
the y- and s-momenta are identical at the injection point, so that 
the natural decrease of y-momentum with increase in s/d. results in 
a decrease in s-momentum during the initial phase of the injection 
process. As the jet axis becomes parallel to the free-stream direc- 
tion, the y-momentum of the jet approaches zero and the x- and s- 
momenta become synonymous . 

The continual decrease of y-momentum along the trajectory of a 
jet injected normal to a cross flow is illustrated in Fig. h9, where 
experimental data of Kamotani and Greber are presented for several 
injection velocities. As you would expect, the y-momentum is largest 
(at a given s/d location) for the Jet with the highest injection 
velocity. These trends of y-momentum with increase in s/d. and 
VR are also reflected by the theoretical results of the present 
study. Abramovich assumed in his analytical development that the 
component of jet momentum perpendicular to the free-stream direction 
(y-momentum) remains constant along the trajectory. The fallacy of 
this assumption is particularly obvious at large s/d i distances. 

Free-stream Nonuniformities 

One advantage of the present theory is its flexibility for 
investigating parameters which affect the trajectory and flow proper- 
ties of the injected jet. Not the least important of these parameters 
are the free-stream velocity and temperature fields into which the 
jet is injected. Up until now the free-stream velocity and 





temperature have been assumed constant, but the next few figures will 
demonstrate some of the effects resulting from relaxing these 
assumptions . 

For the purpose of this illustration, the jet is assumed to 

inject vertically into a free-stream flow which has a boundary-layer 

type of velocity distribution in the y-direction (Fig. 50(a)). 

Similar to the exercise in Appendix A, a Karman-Pohlhausen velocity 

function is described from the injection surface to the boundary— layer 

edge, which is taken to be 8d^. to correspond to the vertical injection 

tests of the present study. At larger values of y/d. the velocity 

is assumed to be constant having the same value as the V used for the 

00 

uniform free-stream velocity case, shown in Fig. 50(b). Figure 51 shows 
that injection into the nonuniform free-stream velocity field results 

in further penetration by the jet into the cross flow than injection 
into the free stream with the uniform velocity field. Coincident 
with this, the jet velocity decay is essentially unaffected, while the 
jet cross-sectional area and momentum are less at any given distance 
along the trajectory. It was noted that the effect of free-stream 
velocity nonuniformity on jet trajectory diminished with increase in 
injection velocity. 

The linear temperature gradient shown in Fig. 52 is used to 
demonstrate the effect of injecting a heated jet vertically into a 
free-stream flow having constant velocity and nonuniform temperature 
fields. The free-stream temperature at the injection surface (T2°F) 
is equivalent to the value used for the uniform free-stream 
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0 < y/d; < 8 V^/V. = 3/ 80 ( y /d.) - 1/51 20 (y/dj ) ^ 

y/d; > 8 v oo /v ± = 0 o 20 

(a) Nonuniform free-stream velocity 



(b) Uniform free-stream velocity 


Figure 50.- Uniformity of free-stream velocity field; a = 90° * 

8 . = 0 °. 1 
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Figure 51.- Effect of fhee-stream velocity nonuniformity 
theoretical trajectory; a. = 90°„ 8. = 0°. 
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Uniform free-stream temperature 


Figure 52.- Uniformity of free-stream temperature field; a, = 90°, 

B. =0°. 1 



1U5 


temperature case. Injecting a Jet with initial temperatures of J2°F 

and 92°F into a free-stream flow having a uniform temperature results 

in the trajectory and temperature curves shown in Fig. 53. The 

combination of the jet flow having a very large Froude number and 

experiencing a rapid heat loss results in no change in the trajectory. 

This result agrees with the experimental observations made earlier in 

the paper as well as with the data of Kamotani and Greber. As noted 

in the figure, the temperature for the heated jet decreases along the 

trajectory until it reaches the free-stream value (T^/T ). Injection 

of the heated jet into the free-stream flow with a temperature 

gradient results in a trajectory similar to that obtained by injecting 

the heated jet into the uniform temperature field. There is a definite 

difference, however, in the temperature curves resulting for these two 

injection conditions, where the nonuniform T situation results in 

00 

higher jet temperatures because of the larger values of T that the 

OO 

jet flow "sees" as it penetrates into the cross flow. Heat is 
initially lost from this Jet flow until a point is reached on the 
trajectory where a heat gain is experienced. 

Three-Dimensional Trajectory 

Up until now all of the theoretical trajectories that have been 
presented are, by definition, two dimensional, i.e. they lie in a 
single plane. The next few figures are presented to illustrate a 
three-dimensional path that results for a jet injected normal 
(cl = 90°) to the mainstream, and rotated k5° away from the vertical 
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when buoyancy is significant. An injection situation similar to this 
was observed in Pig. 31, where the theoretical trajectories are 
actually two dimensional because the Froude number on the jet flow is 
so large, or in other words, the buoyancy force is small compared to 
the momentum forces. In order to obtain a three-dimensional trajectory, 
particularly near the injection point, the Froude number of the jet 
flow must be small. 

For the purpose of this demonstration, the heat transfer from 
the jet is ignored so that the jet density is assumed constant along 
the trajectory. Solutions for the governing conservation equations 
were obtained for Fr. = « (i. e . p.. = pj and for Fr = 10; the 

projections of the resulting trajectories on the x-y and x-z planes 
are presented in Fig. 54. For the Fr. = » casej the proJections 

on the two planes are equivalent which occurs only when 3. = 45° 
and the trajectory is two dimensional. Allowing p. to be less 
than p m leads to the second set of projections, which shows that 
decreasing the Froude number increases the penetration of the jet 
into the cross flow. The fact that the effect is more pronounced on 
the x-y projection than on the x-z projection is an indication that 
the triad of unit vectors associated with the natural coordinate 
system is tracing out a three-dimensional path. In effect these unit 
vectors are twisting" out of the plane where the two-dimensional 
ccto;ry obtained for* the Fz* = 00 case is located « 

Further evidence that these numerical results are consistent 
can be obtained by examining the dependent variables of the governing 





1^9 

equations (u and w) and their gradients. This is done in Pig. 55 
where the injection conditions are identical to those used in Fig. 5^ s 
and the parameters in question are plotted as functions of distance 
along the trajectory. The solid curves represent the case where 
buoyancy force is zero (F^ = 00 ) and are typical of the trends that 
result for injections with 0° £ ol £ 90° and 0° £ £ 90°. As 

noted, as increase in s/d^ results in a continual increase in u 
and decrease in w, u and w approaching 1.0 and 0.0, respectively, 
for large s/d^. This, of course, is coincident with the jet velocity 
vector becoming parallel to the free-stream velocity vector. In 
conjunction with these trends for u and w, du/ds and dw/ds are, 
respectively, positive and negative valued. 

The dashed curves indicate what happens to the trajectory para- 
meters when a sizeable buoyancy force is considered. The fact that 
the values of u and w which result when Fr^ = 10 are, respectively, 
smaller and larger than the corresponding values when Fr^ = 00 is 
indicative of the increase in penetration experienced by the jet when 
the buoyancy force is added. An interesting aspect of this situation 
is seen in the variation of w with s/d^, where w reaches a 
minimum value at some point on the trajectory and then begins to 
increase. This is the same as saying that B achieves a maximum value 

and then begins to decrease. The trend of w with s/d is 

i 

reflected in the dw/ds curve, where dw/ds changes from its usual 
negative sign to a positive value. This result is important because 
it signifies that buoyancy force is becoming dominant over the other 
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forces acting on the jet flow. It is intuitively obvious that for 

the situation where buoyancy is the driving force, the jet path will 

tend toward the vertical so that w will be increasing with increase 

in s/d^ and dw/ds will be greater than zero. This circumstance 

can be shown explicitly by letting = u = du/ds = 0 in the n- 

momentum equation, Eq. (C-9) s presented in Appendix C and by solving 

the resulting expression for dw/ds. 

It is not possible to substantiate the trends discussed here with 

results from other studies , since no experimental data exists for a 

jet with a three-dimensional trajectory. In the theoretical develop- 
29 

ment of Hirst , it was stated that the governing equations could 
apply to a jet following a three-dimensional path, however no calcula- 
tions of a three-dimensional trajectory were presented to support 


that claim. 



XI. CONCLUSIONS 


An investigation has been conducted to determine the important 
parameters that affect the trajectory and flow properties of a heated, 
submerged effluent discharging into a moving waterway . Efforts to 
experimentally and theoretically model the effluent as a Jet injection 
process lead to the following concluding remarks: 

Experimental results for the jet injected laterally, 
vertically, and obliquely into a water channel showed that in- 
creasing jet injection velocity resulted in further penetration 
of the jet into the mainstream. For the range of conditions 
tested, increasing injectant temperature had no discernible 
effect on the jet trajectory for the lateral injection condition, 
decreased the amount of jet penetration for the vertical condi- 
tion, and resulted in a higher trajectory for the oblique injec- 
tion process with the largest injection velocity. Coincident 
with these results, addition of heat to the injectant caused a 
slightly greater rate of spreading of the jet fluid near the 
injection point; this was particularly obvious for the tests 
with the lowest injection velocities. The effects of injectant 
temperature on jet trajectory for the different injection 
conditions were due to a combination of the following factors : 

1) water density is relatively insensitive to temperature change, 

2) the jet fluid experienced a rapid heat loss along the 
trajectory, 3) the injection Froude number was very large implying 
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that the momentum forces dominated the buoyancy forces , 1*) the 
nature of the jet flow was initially laminar for some of the test 
conditions, but tended to become turbulent with increase in 
injectant temperature , and 5) the formation of the counter- 
rotating vortices in the jet flows was different for. the 
respective conditions . 

The theory was developed by using an integral method, which 
accounted for natural fluid mechanisms such as turbulence, 
entrainment 9 buoyancy, and heat transfer, in order to obtain the 
conservation equations governing the jet flow. Solving these 
equations simultaneously yielded predictions of jet trajectory 
and area growth that agreed well with experimental results, and 
thus demonstrated the usefulness of the theory for estimating 
the location and size of the thermal plume with respect to the 
discharge point. 

Unlike previous studies which assumed a specific cross- 
sectional area growth for the jet, the present investigation 
obtained the jet cross-sectional area in the process of solving 
the governing equations. Because of this, the present theory 
provided better estimates for the jet trajectory and allowed a 
prediction for various jet flow properties, such as velocity and 
momenta, to be obtained. Closer agreement with experimental 
jet flow properties could be achieved by improving the an aly tical 
model of the entrainment process. 
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Realistic estimates of temperature in the jet fluid were 
obtained by accounting for heat losses in the jet flow due to 
forced convection and to entrainment of free-stream fluid into 
the jet. Forced convection was seen to be the dominant heat- 
transfer mechanism during the early stages of the jet injection 
process, while the effects of entrainment became dominant as 
the jet penetrated further into the freestream flow. 

The versatility of the theory was demonstrated by observing 
the effects of a jet injected into free-stream flows with either 
a nonuniform velocity field or a nonuniform temperature field. 
Theoretical results were also shown to illustrate a truly three- 
dimensional jet trajectory which was calculated by considering 
the injection Froude number to be small. 
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XIII. APPENDIX A 


Tabulated Test Conditions 

The test conditions associated with the lateral, vertical, and 
oblique injection experiments are presented in Tables III, IV, and V, 
respectively. These test logs contain the values of the parameters 
locating the jet orifice, along with the measured mass flow rates and 
temperatures of the free-stream and injected fluids. Each test is 
assigned a run number to aid in identifying its conditions with the 
photographs presented in the text . The densities p and p were 

OQ 

obtained from the measured values of T^ and T^, respectively, by 
linearly interpolating in a temperature-density table (Ref. 35). 

An average free-stream velocity at the injection station (Vj was 
calculated using 


m = p A V (A-l) 

where A^ is the cross-sectional area of the main flow at the injection 
station. It is noted that has physical meaning only if it is 

considered as the integrated local free-stream velocity over A^, see 
Eq. ( A— 3 ) - The average injection velocity was obtained using, 

*i = P ± A i ^ i (A-2) 

where m is the value listed in the test logs and was obtained from 
the flow meter calibration in Fig.- 17 . 
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TABLE III. TEST CONDITIONS FOP. LATERAL INJECTION PROCESS ; a. = B. ■ 90°, d. - 0.0625 In. 


RUN 

NO. 

a 

inches 

L 

inches 

n. 

1 

Ib^/min 

AT 

i 

o_ 

r 

®oo 

Ib^/min 

T» 

°F 

P ” 

lb f /ft 3 


P i 

lb f /ft 3 

V. 

1 

ft/sec 

0. /p 
1 00 


V v » 





7 

1.25 

27.00 

1.31 

-2.2 

2k2.L 

78.0 

62.266 


62.266 

16.462 


0.9718 

30.095 

30.099 

0.98 

87U8 

_oo 

. 8 



0.79 

-2.2 

2U2.U 

78.0 

62.266 

0.567 

62.266 

9.926 


0.9718 

18.150 

18.153 

0.98 

5276 

-00 

9 



0.23 

-2.0 

21*2.1* 

77-7 

62.268 

0.567 

62.264 

2.S90 


0.9743 

5.283 

5.284 

0.98 

1536 

-387,629 

13 



0.23 

-1.0 

2k2.k 

77.0 

62.256 

0.567 

62.262 

2.890 


0.9870 

5.283 

5.284 

0.97 

1552 

-775,232 


1.90 

26.10 

1.31 

-0.6 

255-6 

75-7 

62.265 

0-577 

62.269 

16.1*61 


0.9921 

28.529 

28.530 

0.98 

8768 


15 



0.79 

-0.1; 

255-6 

75.7 

62.265 

0.577 

62.268 

9.928 

1.00005 

0.9947 

17.206 

17.20T 

0.98 

5276 

-00 

16 



0.23 

-0.4 

255.6 

75.7 

62.265 

0.577 

62.268 

2.890 

1.00005 

0.9947 



0.98 

1536 

_co 

31 • 

1.25 

27.00 

0.79 

6L.5 

2k7.2 

72.0 

62.307 

0.558- 

61.663 

10.061 

0.98613 

1.8958 

18.030 

17.905 

0.53 

98 86 

85,851 

32 



0.23 

5L.2 

2U7.2 

72.0 

62.307 

0.558 

61.612 

2.921 

0.98885 

1.7528 

5.235 

5-205 

0.58 

2623 

9,021 

33* 



0.23 

62.0 

2k7 .2 

71-7 

62.296 

0.558 

61.1*85 

2.927 

0.98698 

1.8647 

5.246 

5.211 

0.54 

2823 

7,746 

34* 



0.23 

63.0 

21*7.2 

71.7 

62.296 

0.558 

61.670 

2.927 

0.98674 

1.8787 

5.246 

5.211 

0.538 

2833 

7,604 

36» 



0.23 

61.3 

2k7.2 

71.7 

62.296 

0.558 

61.450 

2.928 

0.98642 

1.8968 

5.24T 

5-212 

0-53 

2877 

7,427 


•Jet temperatures measured during Runs 33, 34, and 36 using thermocouple probe (see Fig. 11 ) 


TABLE IV. TEST CONDITIONS FOR VERTICAL INJECTION PROCESS; (L » 90°, Sj - 0°, (L = O.O 595 in. 


38 

0 

26.90 

0.21 

-3.2 

240.0 

75.3 

62.268 

0.331 

62.293 

2.911 

1.00040 

0.9575 

8.795 

8.796 

1.02 

ms m 

-251,818 

39 



0.61 

-3.2 

240.0 

75-3 

62.268 

0.331 

62.293 

5.686 

1.00040 

0.9575 

17.172 

17.176 

1.02 



40 



0.69 

-3.2 

240.0 

75-3 

62.268 

0.331 

62.293 

9-565 

1.00040 

0.9575 

28.897 

28.903 

1.02 

■Si* 


45 



0.21 

25.7 

2k0.0 

68.0 

62.326 

0.330 

62.069 

2.922 

0.99591 

1.3779 

8.655 

8.836 

0.80 

1811 

24,786 

46 



0.69 

25-0 

240.0 

69.0 

62.317 

0.330 

62.066 

9.600 

0.99597 

1.3623 

29.091 

29.032 

0.795 

5987 

271,785 

rm 



0.69 

43.7 

222.0 

64.0 

62.352 

0.305 

61.908 

9.625 

0.99283 

1.6828 

31.557 

31.445 

0.69 

6916 

152,950 

PB 



0.61 

44.0 

222.0 

64.0 

62.352 

0.305 

61.906 

5-719 

0.99276 

1.6875 

16.751 

18.683 

0.69 

6109 

53,899 

m 



0.21 

46.0 

222.0 

64.5 

62.349 

0.305 

61.673 

2.931 

0.99237 

1.7132 

9.610 

9.573 

0.67 

2169 

13,318 


TABLE V. TEST CONDITIONS FOR OBLIQUE INJECTION PROCESS; d. = 0.0625 in. 


53 

1.500 


0.69 

9-3 

237.0 

72.0 

62.307 

0.326 

62.228 

8.676 



26.613 

26.597 

0.91 

4965 

706,970 

54 



0.41 

3-6 

237. C 

72.0 

62.307 

0.326 

62.265 

5.152 

0.99933 


15 . 804 

15-798 

0.98 

2736 

4 69,.' ’^9 
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(a. =90°, 

8,-50°) 

0.21 

5-4 

237.0 

70.5 

62.306 

0.326 

62.263 

2.639 

0.99931 

1.0766 

8.095 

8.092 

0.97 

1417 

120,:’6(> 

56 



0.98 

5.4 

237.0 

70.5 

62.306 

0.326 

62.263 

12.313 

0.99931 

1.0766 

37*770 

37.757 

0.97 

6611 

2,618,1 if' 

58 



0.69 

47.6 

237.0 

73.0 

62.286 

0.326 

61.715 

6.748 

0.99083 

1.6548 

26.834 

26 . 711 

0.61 

7669 

98, <A‘> 
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0.41 

44.0 

237.0 


62.286 

0.326 

61.775 

5.193 

0.99180 


15.929 

15.864 

0.63 

4293 

38,HB0 

60 



0.21 

46.4 

237.0 

73.0 

62.286 

0.326 

61.739 

2.662 

0.99122 

1.6356 

8.166 

6.130 

0.62 

2236 

9/i 

61 



0.98 

63. c 

237.0 

73-0 

62.286 

0.326 

61.450 

12.476 

0.98656 

1.8630 

38.270 

38.012 

0.53 

12259 

1 

67 

1.950 

27.50 

0.21 

66.7 

237.0 

64.0 

62.352 

0.305 

61.530 

2.671 

0.98682 

2.0422 

6.757 

8.699 

0.56 

2484 

4, t(«9 


(a, =160°, 

B i =90°) 
















c 9 

1-950. 

26.50 

0.21 

17.1 

237.0 

6s. 0 

62.345 

0.306 

62.205 

2.642 

0. 99775 

1.2631 

8.634 

8.624 

0.90 

1529 


71 

;a 1 =60°. 

Bi-90 0 ) 

0.69 

43.0 

237 .0 

64.0 

62.352 

0.305 

61.916 

5.720 

0. 99301 

1.6719 

28.590 

28.490 

0.695 

6534 

K’H/flil 
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Also shown in the test logs are ratios of temperature, density, 
and velocity between the injectant parameter and its counterpart in 
the free-stream flow. The momentum ratio is an important parameter for 
classifying jet flows and has been used by several researchers'^ 

The square root of the momentum ratio, referred to here as the effective 
velocity ratio (VR) , is seen to be essentially the same as the velocity 
ratio, V^/V^. This is due to the fact that water density is a weak 
function of temperature. 

The other properties of the jet flow at the injection point which 

are listed in the test logs and which were deduced from previous 

information in the tables are viscosity, Reynolds number, and Froude 

number. The viscosity was determined with the values of T from a 

i 

viscosity-temperature plot (see Table A.l in Ref. 36) and as can be 
seen is very temperature sensitive. This dependency on temperature is 
reflected in the values of Re^ and Fr. , which are very susceptible 
to changes in injectant temperature. 

Nonuniform Free-stream Velocity Field 

One of the problems inherent with injection tests performed in 
small scale water channels is the nonuniformity of the free-stream 
velocity field. It is advantageous to investigate the extent of this 
nonuniformity because of its influence on injection characteristics. 

No measurements were made of local free-stream velocities during the 
experimental tests; however, the following exercise was performed 
to shed some analytical light on the subject. 



l6i 


A boundary layer is assumed to form on the channel surfaces as 
illustrated in Fig. 56 (a), where V . is the local free-stream 
velocity. Estimated thicknesses of this boundary layer were previously 
presented in the Tests section of Chapter IX for the flow conditions 
of the current study. The local free-stream velocity can be used in the 
following definition to acquire the average velocity in the channel: 


W. H 




dy dz 


// 


W H 


(A-3) 


dr 


f f 


dy dz 


where If 1 is the cross-sectional area of the water channel flow (A^). 
The idea is to postulate a velocity variation in the boundary layer, 
use the value of yielded by Eq. (A-l), and work backwards through 

Eq. (A-3) to obtain an estimate of the maximum velocity in the channel 
flow. 

Assuming that the local velocity can be divided into regions as 
shown in Fig. 56 (b), and neglecting any shear between the moving water 
surface and the atmosphere , the velocity in the boundary layer regions 
can be approximated by the Karman-Pohlhausen method. This leads to a 
functional form for the local velocity in regions 1 and 2 , where 

= V °o,* (z) ’ ThUS ’ 


= 3 (Z } _ 1 ( z } 3 

V, _ 2 V 2 V 


(A-M 
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where is the velocity at the boundary layer edge and is assumed 

to be constant through region k. The resulting expression for region 3 
is 


JlA = 1 (£) _ 1 ( 2 1) 3 

V 2 V 2 V 

°°.e 


(A-5) 


The local velocity functions shown above are integrated in Eq. (A-3) 
over the respective regions where they are assumed to be acting. 

Taking the proper values for H c ,W^ , and 6 associated with the verti- 
cal and oblique injection tests we find that V = 0.365 ft/sec 

°° s e 

compared to the average velocity value of = 0.330 ft /sec. This 
result suggests that the maximum velocity in the channel is about 10$ 
higher than the average velocity and occupies approximately 75% of the 
flow area. 

The value of for the lateral injection tests would be 

proportionately closer to the value of than was obtained above 

for the vertical and oblique injection tests. This is because the 

estimated boundary layer thickness for the lateral injection tests is 

2/3 of the value used in the calculations above, hence ¥ would 

°%e 

occupy a greater flow area and its value would correspondingly be 
lower. An example of the theoretical effects obtained when the jet 
is discharged into a nonuniform free-stream velocity field of this 
type is presented in Chapter X. 



XIV. APPENDIX B 


Space Curve Information 

This appendix provides additional information about the natural 
coordinate system used in the analytical study , and presents the 
expressions necessary to transform the momentum equations in the text 
to the form used in the numerical solution. 

/ 

Space Curves 

As was mentioned in Chapter VIII, the natural coordinate system 

, , -> 
consists of a moving triad of unit vectors e , e * and e . This 

s n t 

is illustrated in Fig. 57(a) where r is the position vector from the 
origin to a point on the curve and is given by, 

r = x e x + y e y + z e z ( B-l ) 

The derivative of a position vector is shown by Hildebrand^ to be a 
unit vector tangent to the curve and pointing in the direction of 
increasing arc length: 


e 

s 


dr^dx-j . dy ^ dz -* 
ds ds x ds y ds z 


(B-2 ) 


where it is noted that dr/dt is the velocity vector associated with 
a point moving with speed ds/dt along the curve. Thus, 
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is a unit vector, it follows that 


Since e 

s 


+ (&) 
v ds / \ds / 




(B-h) 


The derivative of e with respect to s has a direction 

s 

perpendicular to the curve and is written as , 



d 2 x -*• 

. 2 e > 
ds 


+ & 


ds 


2 y 



( B— 5 ) 


where the length of this vector is the curvature of the curve. Defin- 
ing the radius of curvature (R) as the reciprocal of the curvature, 

-+ , -* 
and e as the unit vector in a direction normal to e , we have 
n s 



1 -*■ 

R ®n 


(B-6) 


so that 



(B-T) 


A moving, rotating triad of mutually orthogonal unit vectors is 

completely described by the addition of the third unit vector, e , 

t 

which by definition is. 
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e t = 


e x e 
s n 


(B-8) 


It is noted that for a plane curve, e and e lie in the plane of 

s n 

the curve, while e is a constant unit vector perpendicular to that 

U 

plane. The trajectory resulting from the vertical injection process 
is an example of this situation. Differentiating Eq. (B-8) leads to, 


d 6 t 

ds 


= -T e 
o n 


(B-9) 


where the scalar T is the torsion of the curve, the negative sign 

o 

implying that T is positive when the vector triad rotates in a 

right-handed sense about e as it progresses along the curve. 

s 

To find de /ds we write e = e, x e and differentiate to get, 
n n t s 



ds 


x e, - 
o t 


1 -* 
R S s 


(E-10) 


Equations (B-6), (B-9), and (B-10) are known as the Frenet-Serret 

formulas. Taking the dot product of e with de /ds leads to an 

w n 

expression for torsion which can be written in determinant form as. 


dx/ ds 

dy/ds 

dz/ds 

d 2 x/ds 2 

d 2 y/ds 2 

d 2 z/ds 

d^x/ds^- 

d^y/ds^ 

d^z/ds 
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We are now in a position where we can find the dot products that 
are required in the momentum equations in Chapter VIII. For example, 
^ e x ’ e n e y * e n ' are needed in the expression for n-momentum, 

Eq. (26). Using equations ( B— 5 ) and (B-6) yields, 


so that we get. 


e 

n 


= R 


d 2 r 

ds 2 


and 




(B-12 ) 


(B-13) 



■ 


ds 2 


(B-lU) 


A similar procedure is used to get the dot products needed for the 
s-momentum expression: 


and 


/■* \ dx 

(e ' e ) = — 

x s ds 


a ■ : > = & 

y s ds 


(B-15 : 


(B-16) 
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The process for obtaining the dot products for the t-momentum 

expression, Eq. (51), is more envolved because of the use of e . The 

"t 

equation for e is found from Eq. (B-10) to be. 


— ~ /r 

T o ds \ ds 1 2 ) 


1 

Ht 


-*■ 

dr 

ds 


(B-17) 


Dotting e^ and 
distributive law 


e , respectively, with this 

J 

for dot products leads to, 


expression and using the 


(t . : } = + r 

x t t ds 2 t 
o ds 


A 

O ds 3 


1 

Rt 


dx 

ds 


(B-18) 


and 



T ds _ c 
0 ds 


1 d£ 

Rt ds 

o 


( B-19 ) 


Direction Cosines 

It was mentioned in Chapter VIII that the procedure for solving 

the governing equations ( 5 ), (2?), (1+3), (52) simultaneously was 

simplified by using the direction cosines of the unit vector e as 

s 

the dependent variables. The angles e makes with the x,y, and z 

s 

axes are shown in Fig. 57(b) and are defined as, 

c °s(a) = cos(B) = ^ cos(y) = ~ (B-20) 
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Eq. (B-U ) represents the auxilliary equation for these direction co- 
sines and is used to express the third direction cosine (dz/ds) in 
terms of the other two. Define 


u 


dx 

ds 


and 



so that Eq. (B-k ) becomes, 


dz _ / 2 2> 1/2 1/2 

__ (1 _ u _ w ) = n 


(B-21) 


(B-22) 


The derivatives of this expression are found to be. 




w 


dw 

ds 




(B-23) 


and 






3/2 



( B-2h ) 


Expressions for R, dR/ds , 

At this point it is desirable to obtain the expressions for R, 
dR/ds, and t q in terms of the direction cosines u and w. 
Equations (B-2l) and (B-23) are substituted into Eq. (B-7) to get the 


expression for R: 
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1 

R 


1/2 


(B-25) 


The relation for dR/ds is obtained by differentiating R, Eq. (B-7), 
and substituting the necessary auxilliary equations from Eq. (B-21) 
to (B-2U): 


2 2 

1 dR _ du d u + dw d v A 

d 3 ds ds 2 ds 2 1/2 

R ds ds n 


(>[■ 


.2 /. \2 "I .2 * 

d w / dw\ A I 

"ds 2 W J " n 3/2 ) 


d 2 u 

U ds 2 VdS ' 


(B-26) 


The torsion of the curve is obtained after some manipulation by 
expanding the determinant in Eq. (B-ll), and by substituting the 
auxilliary equations: 


T o _ / dw du\ ( 1 d 2 u . /du\ 

R 2 = " I" to - " ds) ( n l/2 ds 2 (ds) 


,2 

, d w , 

" ds 2 ^ 


,3/2 


(- 


1/2 


u + n 


1/2 


du\ d 2 w / A 

ds / ds 2 Vn l/2 


w - n 


1/2 dw 
ds 


\ dju 
/ds 2 


(B-27) 


Thus, equations (B-25) to (B-27) represent R, dR/ds, and T q as 
functions of u and w and their derivatives for the general case 
of a three-dimensional curve. It is interesting to note the forms 
these equations take when the trajectory is two dimensional, such as 
for a vertical injection process (y = 90°). For this situation, 
the auxilliary equation (B-22) becomes. 


i 
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2 

u 


+ w 


= 1 


(i .e . n = 0) 


(E-28) 


so that 


A = u 



w 



(B-2 9) 


Using these relations in Eq. (B-25) leads to 


1 _ 1 du 
R w ds 


(B-30) 


Likewise , 


1_ dR - M dfu dw dfw , 

~ 3 ds “ ds , 2 ds . 2 1 

R ds ds 

and T q = 0. Torsion would be expected to equal zero from Eq. (B-9), 
since the trajectory is a plane curve which implies that e is a 

w 

constant vector (parallel to the z-axis). 


0 



XV. APPENDIX C 


Nondimensional Conservation Equat ions 
The purpose of this appendix is to nondimensionalize the governing 
conservation equations and to make several observations concerning 
their numerical solution. In addition s special forms of the governing 
equations are examined and the t-momentum equation is shown to be an 
identity for the vertical injection situation. 


s -Momentum 


The direction cosine expressions from Appendix B are substituted 

. (PAV 2 ) 

into Eq. (13) and the resulting expression is divided by ^ to 


get the nondimensionalized s-momentum equation: 


d. 

1 



where the barred symbols indicate division by d^ , and H = h/d^. 
The expression. 


(pAV) 

E = J - ( C-2 ) 

a. 

1 

was also used in the process of obtaining Eq. (C-l). 

As was mentioned in Solution Procedure in Chapter VIII „ Eq. (C-l) 
is solved at point "j+l" for the jet momentum which is estimated by 


IT 3 




17b 


using a backward finite difference to approximate the rate of change cf 

2 

(pAV ) with s. This gives 


d( PAV 2 ) 
ds 


j+l 


( pAV 2 ) 1+I - (PAY 2 ), 
As 


( C— 3 ) 


which can be substituted into Eq. (C-l)„ The momentum ( pAV 2 ) (pAV 2 ) ^ 
is calculated by using values for the various parameters in Eq. (C-l) 
evaluated at (j) as a first estimate s and then using (j+l) values after 
the first iteration. The Jet momentum is then used in conjunction with 

(pAV) to provide values for 9 and hence A, , . 

J+ 1 j+l j+l 

At this point in the numerical solution the heat loss from the jet 
control volume is accounted for. It should be recalled [see Eq. (53)] 
that the mass of fluid in the control volume must be specified before 
and after entrainment takes place in order to calculate the effect of 
entrainment on jet temperature. The total mass in the control volume 
after entrainment is given by. 


m 


'j+l 


m. + m 

J e 


(C-U) 


where the entrained mass (m ) is found to be, 

e 


m 

e 


(As) 2 E 
V 

a 


(C-5) 


V is the velocity of entrainment and is taken to be an average 

0 , 

velocity in the control volume; Eq. (6) is used in conjunction with 
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the E* function presented in Fig. 3 to specify E. 
n-Momentum 

The expressions for the direction cosines (Appendix B) are 

substituted into Eq. (27) and the resulting equation is divided by 
2 

pAV . The ratio of mass flow in the jet to the initial jet mass flow 
is written as. 


K = 


pAV 

(pAV). 


(C-6) 


and is used to get the nondimens i on al n-momentum equation: 


i = G R ^ + G R 2 (~) + G R % 

R ds Vds^ ds 


where R = R/d. 

x 


and 



( C-7 ) 


In order to put Eq. (C-7) into the form used in the numerical 
solution, the expression for R must be considered [Eq. ( B— 2 5 ) ] . 
Examining the vertical injection case results in a cubic equation 



176 


in du/ds after equations (B-25), (B-28), (B-29), and(B-30) are 
substituted into Eq. (C-7). It is noted that two roots of the cubic 
equation are zero if the buoyancy term is neglected. This type of 
trivial solution can be avoided by dividing du/ds out of Eq. (B-25) 
to get. 



which is substituted into Eq. (C-7) to get the transformed n-momentum 
equation: 


du 

ds 


dw 

du 


+ - 2 . 


(c-9) 


This equation is used to calculate (du/ds), . by knowing u.„ ,w, , 

J+l j+l 0+1 

(dw/du)j +1 , G 1 , G 2 , and G^. Appropriate values of u are then 
obtained by utilizing: 

(l) a backward finite difference scheme at the first step away 
from the injection point (i.e. s = As): 


>i - u j + (f ) 


( C-10 ) 


ds ' j+1 


u 


j+2 


and (2) a central finite difference scheme to get 
on the remainder of the trajectory (i.e. s > As): 


at points 



ITT 


u = u + 2( As ) (— ) ( C-ll ) 

J J Vds/ J+l 

For the vertical or lateral injection situations , the auxilliary 
relation, Eq. (B-22), can he used to obtain a value of w^ +2 
corresponding to the value of Uj +2 in Eq. (C-ll), while its derivative. 


u 


du 

ds 


+ w 


dw 

ds 


0 


( C-12 ) 


can be used to obtain a value of (dw/ds) , and hence (dw/du) . 

J J ' -l* 

These parameters are used to iterate through the governing equations to 
obtain a solution for the two-dimensional trajectory. For the more 
general injection cases, the t-momentum equation is required to provide 
information on w and its derivatives. 


t-Momentum 

The direction cosine expressions from Appendix B are substituted 
into Eq. (52) ana the resulting relation is divided by d^/q^ to obtain 
the nondimens ional t-momentum equation: 

(C-13) 

0 

where x = x d. and 
o o 1 


F X — + F X + F X + F, ( 

1 ° di 2 ° di 2 3 0 U W 


+ F. 


(—) 

\s 2> 


+ F, 


+ f — + F 0 + F — + F x -^ + F x — - + F x = 

T .-2 8 ,-2 *9 - r l0 o .- 11 o —2 12 o 

ds ds ds ds ds ds 
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K A , \ dw 

aT (p ~ - p) ~ 

ii ds 


F = & A R p _ p ) 
2 d.q. ^ p °° p ' 

ii 


F = -£- A ..y. 

3 

d.q.R 

ll 


(p “ P) 


F 


5 




F 


7 





F 


9 



/ u du 
R ds 




We have seen in Eqs . (B-26) and (B-27)that dR/ds and T are 

o 

2 2 2 2 

functions of d w/ds , d u/ds , plus lower order terms. The approach 

2 2 

used here is to define d u/ds by using a finite difference 
approximation after a solution has been obtained from Eq. (C-9). 

The problem then becomes one of getting dR/ds and t q in terms of 

p o 

d w/ds . Thus, from Eq. (B-26) we have. 



and from Eq. (B-27)» 


T 

o 



( c-i4 ) 


( C-15 ) 



179 


where the B coefficients are functions of u, w, du/ds, dw/ds , and 
2 2 

d u/ds . When these relations are substituted into the t-momentum 

P P 

equation, Eq. (C-13), a quadratic expression in d w/ds^ results*. 


B 5 (ft) + B 6 (ft ) 


+ B ? = 0 


(c-l6) 


The coefficients are calculated from information obtained from the 

n-momentum equation, and a solution is acquired at point "j+l" by 
2 2 

approximating d w/ds with a central finite difference. This results 

in a quadratic in w j + 2» which I s solved to provide an update on 

W j+2’ ( dw / ds )j +1 ’ and -(dv/du) . 

It is desirable to demonstrate that the t-momentum equation is an 

identity for the vertical injection situation where the trajectory is 

two dimensional and is confined to the x-y plane. For this case, e 

is parallel to the z-axis which implies that the dot products of this 

vector with e and e in Eq. (51; are zero, and hence the t-momentum 
x y 

equation is an identity. This can be shown in a more rigorous fashion 

• t ^ \ /“K -+ . 

by considering the expressions for (e • e ; and (e ♦ e ), Eqs. 

x t y t 

(B-l8) and (B-19). Substituting for direction cosines into these two 
equations , we can write , 


d 2 u 

ds 2 


dR du u ? 

+ —— + — = o 

ds ds R 


(C-17) 


and 


2 

■o d w dR dw w ? A 

R 72 + i7di* + R = ° 
as 


(C-18) 


c-3 
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where we want to prove that these relations are equal to zero for the 
vertical injection case. To do this it is necessary to use the 
auxilliary expressions that are presented in Eq. (B-22), with n = 0, 
and in Eq. ( C-12 ) . It is noted that Eq. (C-17) can he written as. 


d_ 

ds 



(C-19) 


The expression for R for a vertical injection case was shown in 
Appendix B, Eq. (B-30), and is substituted into the above equation 
to get , 


du , dw ? _ 
u — + w -— = 0 
ds ds 


(C-20 ) 


which is definitely zero from Eq. (C-12). 

The procedure is similar for Eq. (C-l8) where we can write. 



(C-21) 


Eq. (C-12) is solved for dw/ds which, along with the relationship 
for R, is substituted into the above equation to obtain. 


du du 
ds ds 


0 


(C-22) 


Thus , the dot products ( e x 


e^) and (e y 


e^ ) are zero for the 
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vertical injection situation, and hence the t-momentum equation, 

Eq. (5l), is an identity. This fact was used as a check of the 
computer results to see that the numerical output from the t-momentum 


equation was correct. 
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